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We study the production of a single resonance at the LHC and its decay into a pair of Z bosons. 
We demonstrate how full reconstruction of the final states allows us to determine the spin and 
parity of the resonance and restricts its coupling to vector gauge bosons. Full angular analysis is 
illustrated with the simulation of the production and decay chain including all spin correlations and 
the most general couplings of spin-zero, -one, and -two resonances to Standard Model matter and 
gauge fields. We note implications for analysis of a resonance decaying to other final states. 
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I. INTRODUCTION 



■ Physics beyond the Standard Model (SM), to be probed at the LHC PHI], will manifest itself through observations 
of new particles. Such observations are instrumental for establishing the existence of New Physics, though more effort 
is required to understand these observations in detail. It will be crucial to determine the quantum numbers of the 
new particles, their masses, and their couplings to SM fields as accurately as possible. 

Measuring masses, coupling constants, and quantum numbers at a hadron collider is difficult, though many tech- 
niques for doing so were put forward recently. Some of those techniques evolved remarkably over time. For example, 
top quark mass determinations at the Tevatron Q started out from measurements of the tt production cross-section 
and establishing the value of m t which fits the cross-section best. A more recent technique - "the matrix element 
method" - performs a likelihood fit on an event-by-event basis. Since more information about the event is used, 
more efficient separation of signal and background is accomplished and a higher accuracy of the top quark mass 
measurement is achieved. 

The idea that matrix elements, or multivariate per-event likelihoods, can guide us in maximizing the amount of 
information that can be extracted from a given event is appealing; but, to the best of our knowledge, it has not been 
pT} , widely used in hadron collider physics beyond top mass determinations. On the other hand, these techniques are not 
' new to experimental analyses since they were used in many B-physics measurements [5|. 

The goal of this paper is to apply the multivariate likelihood method to the determination of a spin of a resonance, 
produced in hadron collisions. Many extensions of the SM postulate the existence of (elementary) particles of different 
spins that can be single-produced at the LHC. Once produced, these resonances decay into SM particles whose angular 
t-H . distributions contain information about couplings, spins, and other quantum numbers of their parents. 

Spin determination is often discussed in the context of a particular angular distribution; the challenge is to find 
• rH , a distribution that exhibits maximal sensitivity to the spin of a resonance. Single-observable distributions may be 
viable spin-analyzers but, as we illustrate with some examples in this paper, loss of statistical power and certain 
information is inevitable. Construction of the likelihood of the hypothesis that a given event, with its complete 
kinematic dependence, comes from the production and decay of a resonance with a particular spin is the most efficient 
way to analyze the events. Testing this approach in a realistic hadron collider setting is what we would like to do in 
this paper. 

It is best to pursue this program in a situation where the final state is reconstructed fully and accurately. For 
this reason we exclude the final states with missing energy and jets and examine the pure leptonic final states. It 
follows that we can either consider direct decays of resonances to a lepton pair or we can look at the decays of such 
resonances into neutral gauge bosons that subsequently decay into leptons. There are three reasons for us to choose 
the second option. First, more information can be extracted from a fully-reconstructed four-body final state fr3l-fl2l| 
than from a two-body final state; second, direct decays to l + l~ are studied well in the literature 1 |L3l - fl8| : and, third, 
it is reasonable to assume that the decay of a single-produced resonance to Z bosons is sizable, if not altogether 
dominant. Recall that this happens with the SM Higgs boson if its mass exceeds 2m z [l9|. It may also occur in well- 
motivated scenarios of Beyond the Standard Model (BSM) physics. For example, in the extra-dimensional model (20j 



In the appendix we present angular distributions for X — > £+/ that generalize results in the literature. 



FIG. 1: Illustration of an exotic X particle production and decay in pp collision gg or qq — >• X — > ZZ — } 4i . Six angles fully 
characterize orientation of the decay chain: 6* and $* of the first Z boson in the X rest frame, two azimuthal angles <E> and $1 
between the three planes defined in the X rest frame, and two Z-boson helicity angles 9\ and 82 defined in the corresponding 
Z rest frames. The offset of angle $* is arbitrarily defined and therefore this angle is not shown. 

discussed in Refs. [2H - I23I ] KK graviton decays into pairs of gauge bosons are enhanced relative to direct decays into 
leptons. Similar situations may occur in "hidden- valley" -type models (24[. An example of a "heavy photon" is given 
in Ref. H]. 

Motivated by this, we consider the production of a resonance X at the LHC in gluon-gluon and quark-antiquark 
partonic collisions, with the subsequent decay of X into two Z bosons which, in turn, decay leptonically. In Fig. [TJ 
we show the decay chain X — > ZZ — > e + e~ . However, our analysis is equally applicable to any combination of 
decays Z — > e + e~ or It may also be applicable to Z decays into r leptons since r's from Z decays will often be 

highly boosted and their decay products collimated. We study how the spin and parity of X, as well as information 
on its production and decay mechanisms, can be extracted from angular distributions of four leptons in the final state. 

There are a few things that need to be noted. First, we obviously assume that the resonance production and 
its decays into four leptons are observed. Note that, because of a relatively small branching fraction for leptonic Z 
decays, this assumption implies a fairly large production cross-section for pp — » X and a fairly large branching fraction 
for the decay X — > ZZ. As we already mentioned, there are well-motivated scenarios of BSM physics where those 
requirements are satisfied. 

Second, having no bias towards any particular model of BSM physics, we consider the most general couplings of the 
particle X to relevant SM fields. This approach has to be contrasted with typical studies of e.g. spin-two particles 
at hadron colliders where such an exotic particle is often identified with a massive graviton that couples to SM fields 
through the energy-momentum tensor. We will refer to this case as the "minimal coupling" of the spin-two particle 
to SM fields. 

The minimal coupling scenarios are well-motivated within particular models of New Physics, but they are not 
sufficiently general. For example, such a minimal coupling may restrict partial waves that contribute to the production 
and decay of a spin-two particle. Removing such restriction opens an interesting possibility to understand the couplings 
of a particle X to SM fields by means of partial wave analyses, and we would like to set a stage for doing that in this 
paper. To pursue this idea in detail, the most general parameterization of the X coupling to SM fields is required. 
Such parameterizations are known for spin-zero, spin-one, and spin-two particles interacting with the SM gauge 
bosons M, [al and we use these parameterizations in this paper. We also note that the model recently discussed in 
Refs. [2l| - |23l ] requires couplings beyond the minimal case in order to produce longitudinal polarization dominance. 

Third, we note that while we concentrate on the decay X —> ZZ — s- l^l^l^l^, the technique discussed in this 
paper is more general and can, in principle, be applied to final states with jets and/or missing energy by studying 
such processes as X — > ZZ — » X — > W + W~ — > l + vjj, etc. In contrast with pure leptonic final states, 

higher statistics, larger backgrounds, and a worse angular resolution must be expected once final states with jets and 



3 



missing energy are included. We plan to perform detailed studies of these, more complicated final states, in the future. 
However, we note that many results in this paper are applicable to these final states as well. 

The remainder of the paper is organized as follows. In Section [ill we describe the parameterization of production 
and decay amplitudes that is employed in our analyses. In Section ITO1 we calculate helicity amplitudes for the decay of 
a resonance into a pair of gauge bosons or into a fermion-antifcrmion pair; helicity amplitudes for resonance production 
are obtained by crossing. In Section ITVl angular distributions for pp — > X — > ZZ — > /1/1/2/2 for resonances with 
spins zero, one, and two are presented. This is followed by detailed Monte Carlo simulation which includes all spin 
correlations and main experimental effects and which is shown in SectionfVl Analysis using the multivariate maximum 
likelihood technique is applied to several key scenarios to illustrate separation power of different helicity amplitudes 
for all spin hypotheses and in both production and decay, as discussed in Section I VII For completeness, angular 
distributions, including distributions for other decay channels, are given in the appendix. 

II. INTERACTIONS OF AN EXOTIC PARTICLE WITH STANDARD MODEL FIELDS 

In this section, the interaction of a color- and charge-neutral exotic particle X with two spin-one bosons V (such 
as gluons, photons, Z , or W bosons) or a fermion-antifermion pair (such as leptons or quarks) is summarized. The 
spin of X can be zero, one, or two. We construct the most general amplitudes consistent with Lorentz invariance and 
Bose-symmetry, as well as gauge-invariance with respect to unbroken subgroups of SU(3) x SU(2)l x U(1)r of the 
SM. 

The four-momentum of the particle X is denoted by q and the four-momenta of the gauge bosons or fermions 
by qi t 2- The polarization vectors of gauge bosons are denoted by £1.2; we assume them to be transverse q^i = 0. 
Fermion wave functions are conventional Dirac spinors. We employ the field strength tensor of a gauge boson with 
momentum qi and polarization vector ej as /W^" = e?q% — ej'g.f , and the conjugate field strength tensor as = 
1/2 £fj,va(}f^ ,a P — t^vapt-i <zf ■ We use q = q\ — q2 to denote the particular combination of the momenta of the two 
final state particles. 

A. Spin-zero X and two gauge bosons 

The invariant amplitude that describes the interaction between a spin-zero particle X of arbitrary parity and two 
spin-one gauge bosons reads 

A(X H- VV) = v- 1 (>) m ^ +^)/;^ . (l) 

In Eq. ([1]), /* denotes the complex conjugate field strength tensor, v is the SM vacuum expectation value of the Higgs 
field, and A is the mass scale associated with BSM physics. The "couplings" g[°\ 4 are invariant form-factors; the 
upper index reflects the X spin. Since we consider on-shell decays of the particle X to two on-shell gauge-bosons, 
can be thought of as effective dimensionless coupling constants which can, in general, be complex. 

We note that, as written, Eq. {T]) does not use the minimal set of independent variables since it uses both, the 
field strength tensors and polarization vectors for gauge bosons in the final state. However, we write Eq. ([1} in that 
particular way because it can be applied to X decays into both massive and massless gauge bosons and because it 
has the simplest possible connection to SM couplings at tree level. Indeed, if we identify X with the Higgs boson of 
the SM, the proper tree-level amplitude for H — » ZZ is obtained by setting g^ x — and g[ ^ = 2i. To describe the 
coupling of the spin-zero particle to massless gauge bosons (gluons or photons), we simply set 2 m v = in Eq. |T]). 
Clearly, the coefficients gf^ 99 for interaction with gluons, for example, do not need to be equal to the coefficients for 

interaction with the Z bosons gf^' 22 or the photons g^ nl . In fact, Eq. (JTJ) is sufficiently general to accommodate all 
radiative corrections to Higgs interactions with gauge bosons, massive or massless, in the SM, including CP-violating 
form factors that appear at the three loop level 

In spite of the fact that there are four form-factors required to describe the interaction of the spin-zero boson with 
two massive or massless spin-one bosons, there are only three independent structures in the scattering amplitude. To 



2 For X coupling to two gluons, a trivial color factor needs to be introduced in Eq. ( TTl . 
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see this, we rewrite Eq. ([T]) through polarization vectors 

A(X -t VV) = v^e^e™ (axg^ v m 2 x + a 2 q^q v + a 3 e^ af 3 ) i ( 2 ) 
and find the coefficients 01,2,3 to be 

ai=5 (0)l4 + 5 (0)^l + (0) K ^ j a2 = _ 2g (Q)_ 9 (Q) K , a 3 = -2 ff (°). (3) 

m x m x m x 

We have defined the parameters s = qiq 2 = (m 2 x — 2m 2 ,)/ 2 and k = s/A 2 . The amplitude for X decay into two 
massless gauge bosons is obtained from Eqs. @ and ^ by setting m v to zero. 



B. Spin-one X and two gauge bosons 



We consider the case when the exotic particle X has spin one and arbitrary parity. As a consequence of the Landau- 
Yang theorem, the spin-one particle X cannot interact with two massless identical gauge bosons. For this reason, 
a spin-one color-singlet particle cannot be produced in gluon fusion, or decay to two photons. The phenomenology 
of spin-one decays into two Z bosons was recently discussed in Ref. Following that reference, we consider the 
amplitude for the decay to two identical massive gauge bosons X — > ZZ. This amplitude depends on two independent 
form factors 

A(X ZZ) = Mq)(e* 2 e x ) + {e* 2 q)(e\e x )\ + g^e^^e^f. (4) 

Similar to the spin-zero case, g^ and g^ are dimensionless effective coupling constants. We note that these coupling 
constants are, in general, complex with absorptive parts that may arise from quantum loop effects. This possibility 
was not considered in Ref. 9] where the case of zero complex phase difference between the two coupling constants 
was studied. In the case when X has positive parity (J p = 1 + ), the first term violates and the second term conserves 
parity. Alternatively, the two terms correspond to parity-conserving and parity- violating interactions of the l - 
particle, respectively. 



C. Spin- two X and two gauge bosons 

We turn to the spin-two case and construct the most general amplitude for the decay of a spin-two particle X into 
two identical vector gauge bosons. The X wave function is given by a symmetric traceless tensor t^ u , transverse to its 
momentum t^ u q v = 0. Since we would like to apply the formula for the amplitude to describe interactions of X with 
massive and massless gauge bosons, we consider the possible dependence of the amplitude on both the field strength 
tensor and the polarization vectors 



A(X -> VV) = A' 1 [2 5 | 2) t llv f* 1 ' tux f* 2 ' va + 2gi 2 H fll ,^-f* 1 ^ a f* 2 ^ 



A 2 



A 2 



+9. 



(2) 1^ 



A 2 tfiuf 



*l,af) f*(2) 



„( 2 )j- fja 
yio L pa.H 



rq p q°(er(qe* 2 ) + e?(qel)) 



(5) 



(2) 

As in the spin-zero and spin-one cases, g\ t ,, 1B are dimensionless effective coupling constants which are, in general, 

(2) 

complex numbers. They arc different for different gauge bosons V. The first seven constants gl, .,7 correspond to the 
J p = 2 + particle parity-conserving interaction, while the last three terms with <7gf 9il0 correspond to its parity- violating 
interaction. Alternatively, they correspond to parity-violating and parity-conserving interactions of the 2~ particle, 
respectively. 
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We can now write the amplitude through polarization vectors 



+q2fj,qit va ) + c s t Q/ 8^— l-e^po-gfgj + cst af3 qptnv ap q p + Clt 2 ^ {^paq p q a qv + t avp oq p q a qp) 

Tib y ?TZ 
(2) 

The coefficients ci_7 can be expressed through #1,. ,10 

Cl = 2. g p) + 2 g ( 2 )^l + ^) 2 + 2 5 ( 2 )^, 

,( 2 ) 



C2 



C3 



Si 



' (2) 
52 



s 

(2) 



-,(2) 



(6) 



c 4 



-.91 



s 

ml 



ce = 



C7 = gi 'K- 



(7) 



s s 
To describe production of the particle X in hadron collisions, we need to know the X's coupling to gluons. The cor- 
responding amplitude can be obtained from the case A(X — > VV) that we just considered by crossing transformation 

(2) 

and setting m v = 0, gl — 0. Also, because e\q2 = e 2 qi = in the massless case, we find that terms proportional to 
C3 and C4 do not contribute when an analog of Eq. (j6|) is written for massless gauge bosons. 



D. X and two fermions 



For completeness, we also give here the general couplings of the particle X to two fermions. We denote fermion 
masses as m q . We assume that the chiral symmetry is exact in the limit when fermion masses vanish. We obtain the 
following amplitudes 

A(X J=0 -> qq) = ^u qi + p^ l5 ) v q2 , (8) 
A(X J=1 -+ qq) = e%, ( 7 „ (p^ + p^) + ^ + P<Ss)) «», (9) 

A(X J=2 -+ qq) = LlTu* fa (p^ + pfhs) + ! ^ + Pi 2) 7 5 )) tfc, (10) 

where m g is the fermion mass and u and v are the Dirac spinors. It follows that, in the case when fermions are 
massless, the minimal couplings are also the most general ones and no new structures appear. 



III. HELICITY AMPLITUDES 



We are now in position to compute helicity amplitudes for the production and decay processes. Helicity amplitudes 
are important because, as we will see in the following discussion, those amplitudes parameterize angular distributions 
and, hence, can be directly extracted from data. By knowing how these amplitudes are expressed through effective 
couplings introduced in the previous section, we can constrain those couplings through measurements of angular 
distributions. 

To compute the helicity amplitudes A\ 1 \ 2 for the decay X — ¥ VV, we calculate amplitudes presented in the previous 
section for polarization vectors that correspond to Ai, As- We begin with the description of the polarization vectors 
that we use in the analysis. Consider the decay X —> VV in the rest frame of X. The momenta of the two V's are 
parameterized as qi^ — (m x /2, 0, 0, ±/3m x /2), where (3 — (1 — 4m 2 , /to 2 ) 1 / 2 is the velocity of gauge bosons in the X 
rest frame. The polarization vectors for the two Z-bosons read 



< 2 (0) = m- 1 (±/3m x /2, 0, 0, m x /2) , e^±) = e p ( T ) = -L(0, Th -*, 0). 



(11) 
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The polarization vectors of the particle X are defined as follows. For the spin-one boson, Jx = 1, we use 

ex(0) = (0,0,0,1), ejf (±) = ^=(0,Tl,-i,0). (12) 
For the spin-two boson, Jx = 2, the polarization vectors read 

t^(±2) = e x (±)e x (±) , f"(±l) = ^ I4(±)^(0) + e£(0)e£(±)] , 



i""(0) = — [e£(+)e£(-) + 4(-)e£(+)] + J-e^(0)e^(0) . (13) 



It is straightforward to establish general properties of the helicity amplitudes A\ 1 \ 2 . In general there are nine 
complex amplitudes A\ ± \ 2 , since conservation of the angular momentum component along the decay axis fixes the 
spin projection of the X particle to Xx = X\ — A 2 . Because of this identity, we do not reference Xx in the notation 
for helicity amplitudes. Moreover, in the case of two identical vector bosons, such as ZZ, gg, or 77, the number of 
independent amplitudes is reduced from nine to six due to the following identity (2(| [27j 

A^x 2 = (-l) J A X2Xl , (14) 

where J is the spin of the X particle. If parity is conserved, further constraints apply [26|, H3| 

A Xl x 2 = r) P (-l) J A- M -X 2 , (15) 

where r\ P is the parity of the X particle. We note that Eqs. (fl4l [T5)) depend on phase conventions for the polarization 
vectors in Eqs. (fTTl[T2|) . 

For a spin-zero X particle, only Ai — A2 = values are possible. Therefore only ^4++, A , and Aqq contribute. 

For a parity-even scalar with J p = + , such as a SM Higgs, one has A ++ = A and for parity-odd pseudo-scalar 

with J p = 0~ one has ^4++ = — A and Aqq = 0. For a spin-one X particle, Bose symmetry prohibits ^4++, A , 

and Aqq amplitudes, as also evident from Eq. (|14j) . Therefore we are left with only two independent contributions 
A+q = — Aq+ and A-q = — Aq-. Furthermore, in the vector case J p = 1~, one has A+q = A-q = — Aq+ = — Aq~, 
and for an axial vector J p = 1 + the amplitudes are related as A + q = — A-q = — Aq + = Aq- in the case of parity- 
conserving interactions. For a spin-two X particle, there are generally six independent contributions j4oo ; A ++ , 

A , A^ = A |_, A + q = Aq + , and A-q = Aq-. With parity conservation, there are additional constraints for 

the case J p = 2 + : A ++ = A and A +0 = A 0+ = A-q = A -; and for the case J p = 2~: A ++ = —A and 

A + q = Aq + = —^4-o = — Aq-, while Aqq = and A^ = A |_ = 0. We now use explicit expressions for the 

amplitudes constructed in the previous section to illustrate these assertions and compute the independent helicity 
amplitudes explicitly. 



A. Helicity amplitudes for spin-zero decay 

We use Eq. @ to compute the helicity amplitudes and obtain for the decay to two massive vector bosons 

Aoo = -£%^ i(1+/32)+a2/32 )' 

A++ = — { ai + —) ' 
v \ 2 

As expected, these amplitudes satisfy Eq. (fT5j) for parity-even (ai and 02) and parity-odd (03) contributions separately. 

Helicity amplitudes for the decay of the X particle into two massless gauge bosons can be obtained from Eq. (TTB]) 
by disregarding Aqq, since massless gauge bosons cannot be longitudinally polarized, and by setting ft = 1 in A ++ 

and A The relationship between coefficients a\, a 3 and the fundamental couplings gf' , in this case, can be read 

off from Eq. ([3]), where the mass of the vector boson m v should be set to zero. 

It is interesting to point out some features of Eq. (fT6|) for the case when X is the SM Higgs boson, sufficiently heavy 
to decay into two Z bosons. Then, as follows from Eq. (|16p the longitudinal amplitude Aqq dominates in the limit 
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that m H ^> m z . To see this, recall that at tree level in the SM only the term with a\ in Eq. ([2]) contributes. Then, 
A ++ = A__ = — A00/7, where 7 is the boost between the rest frames of the two Z boson 

7 =^L-l = i±^ (17) 

Note also that additional contributions to the helicity amplitudes are present even in the SM, beyond the tree level. 
For example, both a 2 and 03 may appear from radiative corrections. The natural scale for (Z2, generated radiatively 
in the SM,is O(aEw) ~ 10 -2 , while the SM contribution to 03 appears only at three loops and therefore is tiny 
O(10 -11 ) [6jj. In general, all three coefficients ax, a 2 , and 03 are complex numbers with a priori unknown relative 
complex phases between them. In the context of the SM, the measurement of the radiatively induced contributions 
02,3 is a non-trivial test of the Higgs couplings to gauge bosons at the quantum level. 



B. Helicity amplitudes for spin-one decays 



Using Eq. (Q} we find the following helicity amplitudes 



A +0 = -A n+ = 



0+ 



A_ = -A - 



2m z 
I3m 2 x 
2m z 



(18) 



These amplitudes satisfy Eqs. ([HI ITS]) for parity-even (#i ) and parity-odd {g^) contributions separately. Note that 
the Landau- Yang theorem forbids decays of spin-one particles into a pair of massless identical bosons. This feature is 
apparent from Eq. (fT8|) which shows that one of the vector bosons V in the decay of the spin-one particle X should 
be longitudinally polarized. Since massless vector bosons cannot be polarized longitudinally, the decay of X into a 
pair of massless identical vector bosons cannot occur. 



C. Helicity amplitudes for spin- two decays 

Using explicit parameterization for the decay X — ¥ VV in Eq. ^ , we obtain the helicity amplitudes 



A— 



VGA 



V6A 

A +Q = A 0+ = 

A_ = A _ = 

4 

in 

A>o = 



^ (1 + f3 2 ) + 2c 2 /3 2 + i/3(c 5 /3 2 - 2c 6 ) 



ci 



(1 + /3 2 ) + 2c 2 /3 2 - t/3(c 5 f3 2 - 2c 6 ) 



m v \[2A 
3 



,a/2A 



(19) 



As expected, we find six independent helicity amplitudes for the most general case with two identical massive bosons. 
Note that only two independent combinations of 05^.7 constructed from 58.9,10 enter Eq. (|19p . and therefore only those 
combinations are accessible in a measurement. 

To describe decays of X into two massless bosons, all helicity amplitudes with a longitudinal polarization in Eq. (jT9]) 
should be disregarded and (3 — 1 should be substituted everywhere. The relation between the coefficients Cj and 

couplings g\ in this case is found from Eq. J7|) where m v = should be substituted. We note that a peculiar feature 
emerges as the result of this procedure in the case when the X coupling to massless gauge fields is minimal. The 



minimal coupling corresponds to g{ — 1 and g)t,[ — 0. From Eq. ([7]), it follows that in that case C2 = —c%/i = C4/2 



J2) 



and all other coefficients are zero. Eq. (IT9l) then implies that, in the case of the minimal coupling of the particle X 



8 



to massless gauge bosons, A ++ = ^4„_ = 0. Hence, only projections J z — ±2 of the X spin on the collisions axis are 

allowed in that case. However, for more general couplings, A ++ = A amplitudes do not vanish and, therefore, zero 

projection of the X spin on the collision axes is allowed. 



D. Helicity amplitudes for decays into two fermions 

In this subsection, we write down helicity amplitudes for X coupling to two quarks. We need those helicity 
amplitudes to describe production of the resonance X in qq annihilation. The helicity amplitudes for the spin-zero 
case read 



TO, 
V 

The helicity amplitudes for the spin-one case read 



A^ = ^Lm x (p^+pM . (20) 



A+ 



-2m, (> + ^ (pi 1 ) - p£ 
-2m ? (-^ + M( p (x) +Ml ))) 



A_ + = -v / 2m A .(p( 1 )-/3p( 1 )) . (21) 



The helicity amplitudes for the spin-two case read 



A++ = 


2y/2m q m x (3 


V3A 


A— = 


2^/2m q m x (3 


a/3A 


A+- = 


™ 2 X P / (2) 
A \ Pl 


A-+ = 


m 2 x /3 f (2) 
A \ Pl 



(pp'+^o* -/^)) 
^ 2) +#(pi 2) +^ 2) ) 



(22) 



IV. ANGULAR DISTRIBUTIONS 



Information about the quantum numbers of X and its couplings to the SM fields can be extracted from the angular 
distributions. In general, there are five angles that can be studied, see Fig. [1] The two production angles, 9* and 
are defined relative to the parton collision axis; distributions of those angles depend on the production mechanism. 
The three helicity angles, 8\ 1 82, and $, are sensitive to the structure of the interactions with the resonance decay 
products, but they do not depend on the production mechanism. These angles are illustrated in Fig. Q] for the process 
of the type 2—^4, where the decay chain X ~ > Z1Z2 — > (/1/1X/2/2) is considered and the production mechanism is 
cither gg — > X or qq — > X . 

The exact definition of the five angles shown in Fig. [1] is as follows: 8* is the angle between the parton collision 
axis z and the X — ¥ Z\Zi decay axis z' 7 both defined in the X rest frame; $1 is the angle between the zz' plane 
and the plane of the Z\ — > (/1/1) decay in the X rest frame; 9i is the angle between the direction of the fermion /j 
from Zi — > (fifi) and the direction opposite the X in the Zi rest frame, where index i = 1,2 refers to the first or 
second Z boson; finally, $ is the angle between the decay planes of the two Z systems in the X rest frame. The sixth 
angle $* is the azimuth angle of the z' axis with respect to z. It can be arbitrarily defined and it does not carry 
any information about the process. However, $* can be used for bookkeeping purposes to differentiate between the 
first and second Z in the decay in the case that both have identical decay channels. Discussion of the small possible 
difference between the z axis and the beam collision axis follows below. 

Distributions of the five angles described above allow one to determine the spin of the X boson and measure 
contributions of the different helicity amplitudes m both production and decay. The X boson couplings to 
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the SM fields can be deduced from the helicity amplitudes, cf. Eqs. ((TBI - l22l) . We therefore proceed to the general 
expression for the angular distributions in the process pp — > X — > ViV 2 — > (/i/xX/2/2) and illustrate the general 
formula by considering a number of specific examples. We use the helicity formalism as described in Refs. [lEIlSl; f° r 
recent examples of applications of the helicity formalism, see Refs. @, 0, EH, H^-[3l!l . Our results are consistent with 
those references but are, typically, more general. 

Let us first consider the 1 — > 2 decay process. The helicity amplitudes arc defined through a matrix clement of the 
scattering matrix between states with definite projections of angular momenta on a chosen quantization axis. Hence, 



(n, \xMS\Jm) = \J ^£^-D J m * Xl _ X2 (Q)A XlX2 , (23) 

where O describes the polar and the azimuthal angles of one of the final state particles and D denotes the corresponding 
Wigner function. Viewing the collision process of two partons a and b, a b — > X — > ViV 2 — > (/i/xX/2/2)) as a sequence 
of 1 <-» 2 elementary processes, we describe it by the following formula 



A ( XX 0*0 0' 0"\ ( 2J + 1 ) / (2fl + l)(2£2 + l) 

A ab (p a ,pb;xi,X2,m,X 1 ,X2,ni,H2,T 1 ,T 2 ;il ,0,0,0 ) = — — — W — — (24) 

xD£_ X2!m (n*)B xlX2 x D'!: M _ X2 {n)A XlX2 x Dg^WTfa,!*) x DZ* n _ n (tf')W(n, r 2 ) . 

Here the A XlX2 and B XlX2 amplitudes correspond to the X decay and production processes, respectively; T(/ii,/i 2 ) 
and W(ti,t 2 ) describe the decays of the first and the second X daughter to pairs of fermions; Ai j2 are the helicities 
of the two X daughters; Xifl are the helicities of two gluons or quarks in the initial state; [i\,2 arL d Ti 2 are the 
helicities of fermions in the Z decays; and m denotes the helicity of X. We keep the notation general enough with 
J denoting the spin of the X boson and Sx,2 referring to the spins of its decay products, e.g. the Z bosons. We 
choose the convention where O = (0, 0, 0), which by conservation of angular momentum effectively sets m = Ax — A 2 , 
O* = ($>i,d*,— $x), ^' = (0,01,0) and O" = (<&,0 2 ,— Summing over the helicities of all particles, we obtain the 
differential cross-section for the process a b — > X — > V1V2 — > (/i/x)(/2/2) as 

d^dcosfl'd^dcosMcosfl^ = Afab6{Sab " mx)6{ ab x) L I L A ab (p a ,p b] { X ,X;m,n,r};{Q}) ,(25) 

where s a fc is the partonic center-of-mass energy squared, Y ab and y x are the rapidities of the colliding partons ab and 
the resonance X, and J\T a b is the normalization factor 

Nab = 8mlr x c ab (32^r(2m z r z )z ' (26) 



with /3z = yl - 4m|/m^., c g g = 3, c gg = 8, and Tx,z being the decay widths of the resonance X and the Z 
boson, respectively. The relevant differential cross-section for hadron collisions is obtained by convoluting parton 
cross-sections with parton distribution functions 

rftr p p(fl*,$l,flx,6>2,$) = f , , , 7 , s f , , dgqb (XXPX , ^2?2 , 0* , ^1 , gl , gg , g) 

dcos6*d$ 1 dcos6 1 dcos6 2 d3> ^ J x 12 Ja(Xl) Ib[X2> dY x dcos9*d3> 1 dcos6 1 dcos6 2 d$ y a6 =ii„^ ' 1 ' 

ab A x2 

where px,2 are the momenta of the two protons. Note that the angles cos 9*, $x, cos#i, cos 6*2, are defined in the 
X rest frame and are not affected by the integrations over x±, x 2 , and Y x . Convolution with the parton distribution 
functions results in the rapidity distribution of the X boson and may affect angular distributions of decay products 
on an event- by- event basis due to the detector acceptance, discussed later. However, it does not affect angular 
distributions measured with the ideal detector since in this case, integrations over x\ , x 2 , and Y x in Eq. (|27[) factorize 
from all angular dependences. As the result, we may rewrite Eq. (I27[) in the form of the angular distribution of the 
decay products of the polarized particle X 

1 d(7pp _ 1 dr x ^ 28 j 



} pp. tot 



G?cos#*<i$i<icos#ic£cos# 2 <i<i> Tx d cos 9* d<i>xd cos 0\d cos 9 2 d^ 



Since we are interested in the normalized distributions, many normalization factors, as in Eq. (|26p . drop out. However, 
ratios of amplitudes squared | and ratios of partonic luminosities for different production channels appear in the 
normalized angular distributions and contribute to the degree of the polarization of the particle X. Note that if the 
X particle can only be produced in a single partonic channel, all dependence on the partonic luminosities cancels out. 
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This is relevant for both a spin-zero and spin-one X which are produced in ab — gg and qq collisions, respectively. 
For a spin-two X particle, the two parton channels contribute and their relative partonic luminosities affect the 
relative fraction of spin-projection-one polarization as we note below. Apart from this relative normalization, which 
is incorporated into notation discussed below, we can derive angular distributions using the sum of |A &| 2 over the 
helicity states, as indicated in Eq. (|25[) . for each parton channel ab independently. 

To define easy-to-measure parameters related to the helicity amplitudes, we consider the case of a spin-two reso- 
nance as an example. The decay of a spin-two particle to two vector bosons is characterized by nine complex helicity 
amplitudes A\ 1 \ 2 . When vector bosons in the final state are identical, only six helicity amplitudes remain indepen- 
dent. Six complex amplitudes are parameterized by twelve real numbers. However, since we are interested in the 
normalized angular distribution, two of them, the normalization and the overall phase, are unobservable. Hence, 
angular distributions of the spin-two particle decay to two identical vector bosons can be parameterized by ten real 
parameters. We can choose them to be 

A^H^aJVEI^I 2 , (29) 

fc,i=±,0 

0>lA a = ar g(- 4 A 1 A 2 Moo) , (30) 

with (Ai, A2) = (++), ( ), (H — ), (+0), (0—). The remaining (dependent) amplitude parameters fu and 4>ki can be 

expressed either using Eq. (I14[) or the relationship 

ho = 1 - X) /« « ^00 = . (31) 

k,l=±,0 

For the X decays to massless vector bosons or if its spin is less than two, the number of independent non-vanishing 
parameters is further reduced. Our notation for those (simpler) cases are similar to what is described above. For 

example, the spin-zero case is parameterized with four real parameters for (A l7 A2) = (H — h), ( ), and the spin-one 

case with two real parameters for (Ai,A2) = (+0) after re-defining the phase convention with respect to the Aq^ 
amplitude. 

Equivalent parameters can be defined for the helicity amplitudes B XlX2 describing the production mechanism. 
However, not all parameters would enter the angular distributions. We choose the three real parameters f z o, f z i, and 
fz2 to describe the fraction of spin- two X resonance production with spin z-projections 0, ±1, and ±2, respectively, 
where f z o + f z \ + f Z 2 = 1, leaving only two independent parameters. The fractions f z o or f Z 2 arise from the gluon 

fusion mechanism with amplitudes -B++ and B or and B |_ in Eq. (I19p . respectively, where we now use the 

B X1X2 notation in place of A X1X2 to distinguish the process of production from decay. The fraction / zl originates 

from qq annihilation with amplitudes and B |_ in Eq. (j2"2"|) . Note that because we assume that chirality is a 

good quantum number, only ±1 spin projections on the z axis are possible in the annihilation of massless quarks. 
On the other hand, both J z — ±2 and J z = spin projections of a spin-two X are possible in gluon fusion. The 
latter possibility is often ignored in the literature because, accidentally, this contribution vanishes in the case of the 
minimal coupling of the spin-two particle to massless gauge bosons. Finally we note that, since f z \ and f z o,f Z 2 are 
produced by different partonic initial states, those quantities are proportional to corresponding partonic luminosities, 
in addition to the production helicity amplitudes B XlX2 . Below we discuss simplified versions of the general angular 
distribution which are obtained upon integrating over some of the angles in Eq. (|27[) . 



A. Distributions of production angles 

We consider the gg and qq production of an exotic resonance X and its subsequent decay X — > Z1Z2 — > (/i/i)(/2/a)- 
Angular distributions for other decay channels such as X — > 77, X — > gg, and X — > ff are given in the appendix. 

The production angle 6* is shown in Fig. [1] for an X decaying to two Z bosons, as an example. It is defined as 
the angle between the parton collision axis z and the X decay axis in the X rest frame. Determination of the z axis 
requires care since it may differ somewhat from the beam collision axis which we discuss below. The distribution of the 
angle 9* is the only angular observable that contains information about the spin and parity properties of the particle 
X, unless decay chains of the X daughters are analyzed. The latter is possible in the decay to two massive vector 
bosons ZZ or W + W~ for certain cases, but is challenging or practically impossible in other cases. The production 
angle distributions have been extensively studied in the literature [ 1 3l — 1 1 81 ] - However, for the spin-two resonance, we 
point to several modifications of the standard formulas due to generally ignored J z = projection on the collision 
axis. 

Transverse momentum of the X particle introduces uncertainties in the production angle determination. The 
Collins-Soper frame [32| is designed to minimize the effect of X transverse momentum by placing the z axis half 
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way between axes of two beams in the X rest frame. Note that if X has transverse momentum, the two beams are 
no longer collinear in the X rest frame. The uncertainty introduced in the 9* measurement by the non-vanishing 
transverse momentum of the X particle is expected to be relatively small [13] , at least when compared to the statistical 
uncertainties in early LHC measurements. With a larger number of events becoming available later, one can, on 
an event-by-event basis, find and disregard events where X recoils against hadronic jets with a large transverse 
momentum. 

Note that there is an ambiguity in the direction of the z axis and one cannot distinguish 9* from (it — 9*) when two 
identical particles are involved in either the production or decay of the X particle. In particular, this happens if X 
is produced in gluon fusion or when X decays to two identical bosons, such as ZZ, gg, or 77. In those cases, terms 
with an odd power in cos 9* drop out from the angular distributions and we do not consider them further. Only when 
X is produced in qq annihilation and decays into a particle-antiparticle pair, such as W + W~ or , can one try to 
deduce the cos 9* sign. 

The case of a spin-zero particle is very simple. Since no spin correlations are involved, the production angles 
distributions are flat. However, determination of these angles is still relevant for the analysis of a spin-zero particle 
because it is needed to discriminate against backgrounds or resonances with non-zero spin. Normalized distributions 
of production angles for spin-one and spin-two resonances that decay as X — > ZZ — > (/i/i)(/2/2) can be written as 
follows 



327rdr 



3rdcos6»*d$i 



2(1 + cos 2 9*) - V2/+o(l - 2/ +0 ) sin 2 9* 008(2$! - [0 +o - &>-]) • 



(32) 



k^J^L'L = ( 2 - 2 ^ + - 6 ( 2 - ~ /-a) cos2 r + 3 ( 6 - 10 /*i - cos 4 9* 
51 acos9*d<Pi 

+f+-{(2 + 2f zl - 7f z2 ) + 6(2 - 6/ 2l + f z2 ) cos 2 9* - 5(6 - 10/ 2l - 5/ 22 ) cos 4 9*} 
-2(/ +0 + / _) {(2 - 4/,i - f z2 ) - 6(4 - 7f zl - 3/ 22 ) cos 2 9* + 5(6 - 10/ 2l - 5/ 22 ) cos 4 9*} 
-2v7+o/o-{(/*i - fz2) + (6 - 10/ 2 i - 5/ 22 ) cos 2 } sin 2 9* cos(2$! + +o - 4>o-) 
-~ {(2 - 2f zl - 3/ 22 ) - (6 - 10/ 2l - 5/ 22 ) cos 2 9*} sin 2 9* 



x 

3?ri?i 



4 

x 



VQf++f+- C0S(2$! - ++ + (/>+-) + y/6f— /+-C0S(2$ 1 4 

{(2 - 4/ 2l - f z2 ) - (6 - 10/ 2l - 5/ 22 ) cos 2 9*} cos 9* sin0* 



V3/+o(l - /++ - /— - V+- - 2/+o - 2/o_) cos($ 1 - ^ +0 ) + v /3/ +0 / ++ cos($ 1 + +o - 0++) 
- V / 3/o-(l - /++ - /— - 2/+_ - 2/ +0 - 2/o_ ) cos($ 1 + 0o-) - V3/o-/— cos($! - O - + 0— ) 
■^^{(6 - 6/ 2 i - 9/ 22 ) - (6 - 10/ 2 i - 5/ 22 ) cos 2 9* } cosfl* sin^ 

h- - 0+0 ) - v / 2/+-/o-Cos($i - 0+_ + 0o- 



^2/+_/ +0 cos($i 



(33) 



Note that the right hand side of Eq. (|33[) is written as a linear combination of terms of the form («i +a 2 / 2 i +cK3/ z2 ), 
multiplied by other parameters. It is peculiar that in all such terms, except in the very first, 9*- and $1 -independent 
one, there is a relationship between the a-coefficients ot\ = —0.4 x (a 2 + 0:3). This means that if we choose f z \ = 
f z2 = 0.4 in Eq. (|33[) . the production angle distribution for the spin-two particle becomes flat! To understand this, 
note that those helicity fractions imply production of an unpolarized X-boson in which case the production angle 
distribution must be constant. 

It is interesting to point out in this regard that observation of the flat production angle distribution in a two- 
body decay X — ¥ P\P 2 does not mean that the X's spin is zero. In fact, X can have any spin J but, for the 
flat distribution, it must be produced unpolarized. In turn, this implies that all helicity fractions must be equal 
2/ 2 o — fzi — jzi — ■■■■fzj- In general, each helicity fraction f zm is a sum of many terms 



\9» 



b^X,m\ 



*fafb 1 



(34) 



ab 



where a, b are the two partons whose collision produces the resonance X, g a b^x,m are the couplings of the two partons 
to the resonance X in the helicity state m and f a & are the parton distribution functions. It follows that the equality 
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of all helicity fractions requires an unnatural tuning between the coupling constants and the parton distribution 
functions. For example, a spin-one resonance can only be single-produced in qq collisions and, in general, there are 
two helicity fractions f z o and f z \. Since, as follows from Eq. (|2"Tj) . f z0 is proportional to the quark mass squared, it 
is very unnatural to expect unpolarized production in the J — 1 case. For an X particle of spin two, unpolarized 
production requires tuning of the coupling constants and parton distributions since gluon collisions are responsible 
for the f z o and f Z 2 helicity fractions and qq collisions - for the f z i helicity fraction. It is interesting that for a fixed 
mass X resonance, such tuning can only be argued for a particular energy of the hadron collider; changing the collider 
energy, at least as a matter of principle, will clearly destroy the tuning and, if J > 0, will turn flat cos0* distributions 
into non-flat. While the above discussion shows that unpolarized production of the resonance with non-vanishing spin 
requires a high degree of tuning between coupling constants and parton distribution functions, we emphasize that, as 
a matter of principle, observation of a flat production angle distribution does not immediately imply that the spin 
of the resonance is zero. As we discuss next, the analysis of helicity angle distributions helps in distinguishing the 
different spin scenarios. 



B. Distributions of helicity angles 

It follows from Eq. (|25l) that the most general angular distribution depends on five angles. Such a distribution 
contains information about production of the resonance X and its decay into ZZ — > {fifi){f 2 f2)- In this section, we 
restrict the presentation to angular distributions which are averaged over the production angle 3>\. The most general 
angular distributions are given in the appendix. To obtain distributions differential in the three helicity angles, one 
can easily integrate over the cos 0* in the formulas below. In a similar manner, it is easy to obtain one-dimensional 
projection of any of the five angles that describe the decay. In this subsection, we present those distributions for 
spin-zero, spin-one, and spin-two resonances that decay asl^ ZZ —> {fiji){f2h)- Because Z decays are involved, 
distributions depend on the parameters R1.2 = 2ri j 2/ / (l + r\ 2 ), where 7*1.2 is the ratio of axial to vector couplings 
of the fermions f\.2- Specifically, r — ca/cv — t 3L /(t 3L — 2<7sin 2 9w), where q is the fermion charge and t 3L is its 
weak isospin. In this paper we mostly consider Z decays to charged leptons. For them, q = —1, t 3L = —1/2, and 
R ~ 0.15. Note that if the first Z decays to charged leptons and the second Z decays to quark and antiquark jets, 
one needs to average the value of R2 over contributing quark flavors since R is different for up and down quarks, 
i? up ~ 0.67, -Rdown — 0.94. However, if the quark and antiquark jets cannot be distinguished in the Z decay, this is 
equivalent to R 2 = in the measured angular distributions. 

The angular distribution for a spin-zero resonance is independent of the production angles. It depends on four free 
parameters and reads 

4(1-/++-/—) sin 2 9 X sin 2 6 2 

+ (/++ + /— ) ((1 + COS 2 6x)(l + COS 2 82) + 4i?ii? 2 COS*?! COS0 2 ) 

-2 (/++ - /__) (Ri cos 0i (1 + cos 2 2 ) + R 2 (l + cos 2 Ox) cos0 2 ) 
+4V/++(1 - /++ - /— ) (Ri - cos0i)sin0i(i? 2 - cos0 2 )sin0 2 cos($ + (/)++) 
(1 -/++ -/— ) (Ri + cos0!) sin 0i (-R 2 + cos 2 ) sin 2 cos($ - <£__) 
+2 v //++/__sin 2 0i sin 2 2 cos(2$ + 0++ - ) . (35) 

We point out that non-zero values of Ri are reflected in preferential directions of fermions in Z decays, see e.g. terms 

i?icos0i and i?2COS02 that are present in Eq. (|35|) if parity is violated with /++ 7^ / Equation (|35|) is the most 

general angular distribution of the decay of a spin-zero particle and, as such, generalizes many similar results presented 
in the literature. 

We note, however, that in specific cases the number of parameters in Eq. (|35j) can be reduced. For example, 
considering the tree-level coupling of the SM Higgs boson to two Z bosons, we find the relation between three helicity 
amplitudes to be A ++ = A = — Aqo/j, where 7 is defined in Eq. (fTT|) . This leads to 

/++ = /— = 2— , 0++ = <t>-- = 7T • (36) 

It is easy to account for changes in H — > ZZ — > Al angular distributions caused by the radiative corrections in 
the SM, if of M and af M are known. Similarly, for a parity-conserving interaction of a J p — particle, one has 
/++ = / — = 1/2 and <j) ++ - cj> — = tt. 



1287rdr J=0 
9T(icos0idcos 9 2 d<^ 
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The angular distribution for a spin- one particle decaying to X — > ZZ is determined by two parameters; we choose 
them to be f +0 and [^ +0 — 4>a-] = ar g04+o/A)-)- The angular distribution reads 

o ?T ^ 5 I 2 I dr l^ ZT^ = (l + cos 2 r| x {l-cos 2 1 cos 2 e 2 + (l-4/ +o )(i? 1 cos^ 1 sin 2 2 + i? 2 sin 2 1 cos0 2 ) 

27 T d cos 8* d cos Old cos0 2 a<P I J I 

+ v/8/+ (l - 2/ +0 ) sin0i sm9 2 (R 1 R 2 - cos0i cos0 2 ) cos($ + [0 +o - O -])} . (37) 

This distribution differs from a recent result presented in Ref. [9(. The origin of the difference is the reality of the 
two effective couplings assumed in that reference. Indeed, if the two couplings are real then, independent of 
their actual values, / + o = 1/4, as follows from Eq. (fT5|) . Then, the terms in Eq. (|3"T)l that are linear in Ri do not 
contribute. In the more general case of two complex couplings, the preferential polarization of the Z boson appears 
and gets reflected in the preferential direction of, say, the negatively charged lepton in Z-decays. We also note that 
the angular dependence in Eq. (|37|) is a product of the production angle distribution and the distribution of three 
helicity angles. 

Finally, we consider a spin-two resonance. The normalized angular distribution is 



2048?r dTj =2 



45 Td cos 0* d cos 9i d cos 9 2 d§ 



{(2 - 2f zl + f z2 ) - 6(2 - 4/ 2l - f z2 ) cos 2 6* + 3(6 - 10/ zl - 5f z2 ) cos 4 , 



x {4 (1 - /++ - /__ - 2/+_ - 2/ +0 - 2/ _) sin 2 X sin 2 9 2 

+ (/++ + /—)((! + COS 2 00(1 + COS 2 9 2 ) + ARlR% COS01 COS 9 2 ) 

- 2 (/++ - /__) (R 1 cosfl^l + cos 2 6 2 ) +R 2 (1 + cos 2 0i) cos9 2 ) 

+ - /++ - /— - 2/+_ - 2/ +0 - 2/ _)(i?i - cos 0i ) sin X (R 2 - cos 2 ) sin 9 2 cos($ + 0++) 

+ V/— (1 - /++ - /— - 2/+_ - 2/ +0 - 2/o-)(.Ri + cos0i) sin0i(i? 2 + cos0 2 ) sin0 2 cos($ - 0__) 
+ 2 v // ++ /__sin 2 0i sin 2 2 cos(2$ + 0++ - 0—)} 

+8{(/,i + f z2 ) + 3(2 - 3/ 2i - 2f z2 ) cos 2 0* - (6 - 10/ zl - 5/, 2 ) cos 4 0*} 

X {(/+o + /o-)(l - cos 2 0i cos 2 2 ) - (/+„ - f -)(Ri cos0i sin 2 2 + R 2 sin 2 Q x cos0 2 ) 
+ 2\//+o/o- sin 0i sin0 2 (i?ii? 2 - cos0i cos0 2 ) cos($ + [tp +0 - 0o-])| 

+ {(6 - 2/ 2i - 5/rt) - 6(2 - 2/,i - 3/, 2 ) cos 2 0* + (6 - 10/,i - 5/ 22 ) cos 4 0* } 

x /+_{(1 + cos 2 0i)(l + cos 2 2 ) - 4i?ii? 2 cos0icos0 2 | , (38) 

where the dependent parameter / o is expressed as (1 — / ++ — / — 2/^ — 2f+o — 2/o_). 

There are two ways to obtain an angular distribution in the three helicity angles from Eq. p8l) . One way is to 
integrate over the cos0* in Eq. (13"8"|) ; note that Eq. (|3"8")) is written in such a way that the dependence on cos0* 
in each term is factored out. Upon integrating over cos0*, the dependence on f z i and / z2 must disappear since 
the distribution over helicity angles is independent of the production mechanism. We note that this independence 
suggests an alternative way to obtain the distribution of the helicity angles from Eq. (f3"5]) . Imagine that X is produced 
unpolarized, which amounts to setting f z i — f z2 = 0.4 in Eq. (j3"8")l . In this case coefficients of all cos 0*-dependent 
terms in Eq. (|38p vanish and integration over that angle becomes trivial. 

Even if X is produced unpolarized, it is still possible to separate spin hypotheses by analysis of the helicity 
angles. The opposite situation is also possible when helicity angles provide no separation while productions angle 
distributions are different. In fact, joint analysis of production and helicity angular distributions is the most efficient 

way to separate different hypotheses. However, there is a special case of the helicity fractions /_| = / + o = fa- = 

and production polarization j z \ — f z2 — 0.4 where the angular distributions in Eqs. (|35|) and (|38[) become identical. 

The situation /_| = / + q = /o_ = is realized when only terms with c 2 , C3, and C5 contribute to helicity amplitudes 

in Eq. ([5]), which becomes equivalent to Eq. @. This is possible, though very unnatural, situation where the 

amplitude contributes in production but vanishes in the decay mechanism. However, this possibility cannot be ruled 
out experimentally unless different spin projections in the production arc probed. 

It is a general feature of Eq. (j3"B"|) written for spin-two decays that it includes both Eqs. (|35p and (|3T|) written 
for spin-zero and spin-one, respectively. Note that Eq. (13"8"|) factorizes into the product of a function of cos 0* and a 
function of three helicity angles under the following three conditions. In the case /_| = /+o = fo- = 0, the helicity 
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angular distribution is given by Eq. (l35j) and the production angle distribution appears in Eq. (|33|) . This trivially 
reduces to the case of unpolarized production discussed above. In the other case /o_ + /+o = 1/2, and therefore the 
other fki = 0, the helicity angular distribution is given by Eq. (|37]) and the production angle distribution can be read 
from Eq. (1331) . In this case it happens that, under the condition f z \ = 0.6 and f Z 2 = 0, the function of the production 
angles is proportional to (1 + cos 2 9*) and Eqs. (l37|) and (|38|) become identical. However, this is only a coincidence 
and the full angular distributions which also include angle $i do not match, as for example Eqs. (|32|) and (|33|) arc 

not identical in this case. Finally, the last term in Eq. (|38p corresponds to the non-zero A_| = A |_ amplitudes and 

is unique to the spin-two case, though its helicity angular distribution is very similar to the other transverse term 
which remains in the case /++ = / = 1/2. 



We have written a Monte Carlo program to simulate the production and decay of spin-zero, spin-one, and spin-two 
resonances in hadron-hadron collisions, including all spin correlations, in the processes gg(qq) — > X — > ZZ — > 4L A 
special feature of our program, that distinguishes it from other recent implementations 33], is that it includes the 
most general couplings of the X particle to gluons and fermions in production and to Z bosons in decay, as explained 
in the previous sections. Extension to other final states, including hadronic Z decays, is straightforward. 

The spin-zero resonance is produced in collisions of gluons, and the spin-one resonance is produced in qq collisions. 
For the spin-two resonance, both partonic channels are included in the simulation. Since, in the general case, the 
relative strength of ggX and qqX couplings is not known, the program allows the request that the spin-two X particle 
is produced with fixed relative frequency in qq and gg collisions. We note that both the resonance X and the two 
Z bosons are considered to be on-shell in our program; all the off-shell effects are neglected. The narrow width 
approximation is a reasonable assumption for the Z bosons, but it is not possible to say a priori if this also is a good 
approximation for the resonance X. The applicability of the narrow width approximation for the particle X is the 
assumption built into our program, but it is relatively straightforward to remove this restriction in the future. 

The program can output both weighted and unweighted events, depending on the requested mode of operation. 
Weighted events are typically used for fast calculations of simple one-dimensional distributions and for debugging 
purposes. Unweighted events, on the other hand, are used to interface the results of our program to programs that 
simulate realistic detector response. 

We now turn to the discussion of how backgrounds of the resonance X for four-lepton production in hadron collisions 
are simulated. Such backgrounds arise from Z{ r y)Z( r y), Zbb, ti, W + W~bb, WWZ, tiZ, 4b etc, where we assume that 
6-quarks decay semileptonically. Backgrounds that involve ^-decays into leptons can be controlled by requiring that 
leptons are isolated. Other backgrounds can be strongly suppressed by requiring that the invariant masses of lepton 
pairs are close to the Z-boson mass, that there is no missing energy in the event, and that the four tracks originate 
from the same vertex near the interaction region. As a result, pp — > ZZ — » 41 is the major irreducible background 
that survives all the possible selection requirements [34l. 135} . We emphasize that, while it is possible to understand 
gross features of the backgrounds, it is conceivable that subtle (but relevant) details of the angular distributions due 
to background processes may weakly depend on the exact background composition. Such an exact composition is hard 
to predict theoretically with any degree of confidence but, fortunately, this is not necessary. Indeed, this problem can 
be solved by using sideband analysis. This possibility is incorporated into our analysis discussed in the next section, 
but its detailed study is beyond the scope of the present pap er. Here, we restrict ourselves to pp — > ZZ — > 41 as the 
only background source and simulate it using Madgraph [36j. 




In order to illustrate MC simulation, compare it to the derived angular distributions, provide examples of the 
data analysis techniques, and understand the statistical power of the proposed approaches; we choose seven scenarios 
which cover all spin and parity combinations. They are described in Table |TJ Note that for the 2 + graviton-like 
resonance, we consider two models with different couplings to matter and gauge fields. The distributions of all five 
angles - cos 0* , $1, cos#i, cos 62, and <& - for the seven models in Table Q] and the pp — > ZZ — > 41 background 
process are shown in Fig. [21 These distributions were generated with our program assuming that the resonance mass 
is m x — 250 GeV. Throughout the paper we consider ^/s — 14 TeV proton-proton collisions and use the CTEQ6L1 
parton distribution functions [37], H|| . Note that the distributions of cos Q\ and cos 82 are identical and are combined 
in one plot in each case. Projections of the ideal angular distributions derived in the previous sections agree well with 
simulated distributions. A glance at Fig. [2] suggests that different hypothesis about resonance quantum numbers can 
be efficiently separated if all five angles are analyzed simultaneously. Of course, correlations in the multi-dimensional 
space of all angles are important for full separation power and those correlations cannot be easily illustrated. We fully 
exploit those correlations in the angular analysis discussed in the next section. 

Our MC generation is performed stand-alone but, since unweighted events are produced, it is easy to incorporate it 
into a software framework that includes full detector simulation. This is achievable in the same way as for Madgraph 
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TABLE I: The list of scenarios chosen for the analysis of the production and decay of an exotic X particle with quantum 
numbers J p . For the two 2 + cases, the superscripts m (minimal) and L (longitudinal) distinguish two scenarios, as discussed 
in the last column. When relevant, the relative fraction of gg and qq production is taken to be 1:0 at m x = 250 GeV and 3:1 at 
m x = 1 TeV. The spin-zero X production mechanism does not affect the angular distributions and therefore is not specified. 



scenario (J p ) X — > ZZ decay parameters X production parameters comments 



0+ 


a% ^ in Eq. © 






gg^x 
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X: p 21 , p 22 / in Eq. (0 


"pseudo-tensor" 



interfaced through Pythia 39]. However, to illustrate effects of realistic detector response, we employ a simplified 
technique not attached to any particular experiment. 

Note that, with our choice of the final state, we require measurements of the four-momenta of all charged leptons 
for complete reconstruction of the event kinematics, including boosts to the rest frames of X and Zs, where the 
production and helicity angles are defined. In an experiment like ATLAS or CMS, the four-momentum of the charged 
lepton is reconstructed from its track in the detector and there are two main effects that influence these measurements: 
(1) deviation of the five measured track parameters, transverse momentum (j>t), direction (4>, 9), and distance to 
the origin (d xy , d z ), from their true values, and (2) non- uniform efficiency of particle detection across geometric and 
kinematic parameters. We model the detector response in the following way. First, we assume that, for typical track 
momenta of order 50 — 250 GeV/c considered in our examples, the track parameter resolution can be identified with 
the typical ATLAS or CMS tracker system resolution. Throughout this paper, we model detector resolution using the 
recently reported CMS track resolution parameters [4(| , obtained from analysis of cosmic ray data. For the three track 
parameters {pt, <t>, we apply conservative Gaussian random smearing with an rms Apx = 0.025 x px + 0.0001 x p\ 
(GeV/c), A<fi = A6 = 0.001 (rad), and neglect resolution effects on the track origin. Within this simplified model of 
detector response, we observe that production and helicity angles can be measured with a typical resolution of the 
order of ~ 0.01 rad. This resolution is rather good and does not infringe on our ability to perform angular analysis. 
Note that because production and helicity angles are defined in Lorentz frames which differ from the laboratory frame, 
they are affected by uncertainties in all track parameters and, most importantly, by pt uncertainties. 

The second detector effect is the non-uniform reconstruction efficiency. We model it in a simple way assuming that 
only tracks produced centrally can be measured. We also assume that the detection efficiency does not change within 
the detector acceptance. Hence, the acceptance function is given by the step-function 

0(0*,*1,01,02,*)= II %max|-|li, fl («*»^,«l.«2»^)l), (39) 
i=l,2; q=±l 

where r\i.± — lncot(#i±/2) is the pseudorapidity of a lepton with the charge q = ± that originates from the decay of 
the i'th Z boson, |?7 m ax| = 2.5 is the maximal pseudorapidity that the tracker can reconstruct, and the product runs 
over the four leptons in the X decay. 

Given the simple form of the acceptance function, it is straightforward to understand the effect of detector acceptance 
on angles measured in the laboratory frame, but their effect on the production and helicity angles is far less obvious. 
Indeed, the rapidity of the lepton i is a function of all production and helicity angles; as the result, Eq. ([3"9"| cuts out a 
complicated region in the five-dimensional space of production and helicity angles and makes the efficiency dependent 
on those angles. For example, if X and the Z bosons are produced at rest in the laboratory frame (this, of course, 
requires m x = 2m z ), the polar angles of the four final state leptons #i+,i-.2+.2- that appear in Eq. (|39p can be 
expressed through the five angles characterizing the decay as (with the convention $2 = $ + 

=F cos 8i± — cos 9* cos 6*; + sin 6* sin 9{ cos <F; . (40) 

However, in general, X is boosted in the laboratory frame and the Z bosons are boosted in the X frame. This 
introduces an additional transformation of the four 9i± angles. We illustrate the effect of the lepton rapidity cut on 
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FIG. 2: Distribution of the cosf?* (left), $i (second from the left), cos#i and cos 82 (second from the right), and $ (right) 
generated for m x = 250 GeV with the program discussed in the text (unweighted events shown as points with error bars) and 
projections of the ideal angular distributions given in the text (smooth lines). The four sets of plots from top to bottom show 
the models discussed in Table U for spin-zero 0+ and 0" (top), spin-one 1 + and 1 (second row from top), spin-two 2^, 2j, 
and 2~ (third row from top), and the bottom row shows distributions in background generated with Madgraph (points with 
error bars) and empirical shape (smooth lines). The J + distributions are shown with solid red points and J~ distributions are 
shown with open blue points, while the 2„ and 2j are shown with red circles and green squares, respectively. 



production angles in Fig. [3J where we plot the distributions of 9* and <&i production angles for the spin-zero particle 
X. If these distributions are measured with the "ideal" (Air) detector, the results are flat. Hence, the non-trivial 
shapes of these distributions shown in Fig. [3] are entirely due to an acceptance effect. 

It is evident from Fig. [3J that the acceptance effects are very important in the analysis of data. They have to be 
taken into account explicitly, otherwise the results of the analysis will be biased. This can be easily done in our MC 
simulation program on an event-by-event basis using the acceptance function in Eq. ([39]), where we reject events if 
at least one lepton exceeds the maximal pseudorapidity. It is also possible, but much harder, to incorporate this 
acceptance function into the likelihood function that is discussed in the next section. However, as we explain now, 
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C0S6* <!>., 

FIG. 3: Distribution of the cos 6* (left) and $i (right) for the case of spin-zero resonance production gg — > X — > ZZ. The 
mass of the resonance is m x = 250 GeV (solid points) and 1 TeV (open points). Detector acceptance effects are taken into 
account, see text for details. Lines show empirical parameterization. 

certain simplifications in the treatment of the acceptance function are possible. These simplifications allow us to 
derive a probability distribution that accounts for major acceptance effects. 

The first of the simplifying assumptions is easy to understand. Suppose X is a heavy resonance. Then, Z bosons 
from X decays are strongly boosted and decay into highly collimated "leptonic jets". Therefore, as the X mass 
increases, in the laboratory frame the decays look exceedingly pencil-like, with all the leptons aligned with directions 
of their parent Z bosons. It is then clear that the detection efficiency of the leptons is highly correlated with the 
Z-production angle 0* but not with the other angles. If the X boson has an intermediate mass, the story becomes 
more complex and acceptance effects in other angles should be expected. However, it is reasonable to expect those 
effects to be moderate. We can check these features of the acceptance with the MC simulation. From Fig. [31 it is 
evident that with the increase in the mass of the particle X the acceptance effect in cos 0* gets more pronounced 
while it is indeed reduced in the other angles. 

The second simplifying assumption is based on the observation that distributions of the helicity angles 
(cos 0i, cos 02, are independent of the production mechanism. Hence, there are no acceptance effects in the helicity 
angles provided that cos 0* and $i are random, as in the case of spin-zero particle production. In a more general 
case, acceptance effects on the helicity angles are expected to be rather mild. 

VI. DATA ANALYSIS 

In this section, we illustrate the application of the full angular analysis formalism by considering the production 
of a resonance X and its subsequent decay to two Z bosons. We use the MC simulation discussed in the previous 
section to generate both signal and background events, and apply a multivariate fitting technique discussed later in 
this section to extract as much information about the produced resonances as possible. We point out that, because 
our analysis is general and because it is not based on any particular model of BSM physics, it is not possible to 
accurately predict how many X particles are produced, for a given luminosity. Instead, we find it more reasonable 
to assume a certain number of signal and background events, reconstructed by an LHC experiment. This number of 
events should be sufficiently large to enable meaningful angular analysis, yet it should be small enough, to be relevant 
for the situation soon after the resonance observation. The statistical power of the method can be easily extrapolated 
to a different number of resonance events. 

As we explain below, we perform two different analyses in this paper, for which we choose to have 30 and 150 signal 
events, respectively. Roughly, this corresponds [HI, IH[ to the number of SM Higgs boson events pp — > H — » ZZ — > 41^, 
I = e, (i, detected by one LHC experiment with 5 or 25 fb -1 integrated luminosity, for a Higgs boson mass m H = 250 
GeV, at a proton-proton center-of-mass energy of 14 TeV. Of course, as we pointed out a few times already, the 
actual production rate of exotic particles may be larger or smaller; we choose 30 and 150 events merely for illustration 
purposes. 

We consider two mass points, the low mass m x = 250 GeV and the relatively high mass m x = 1 TeV. The 
background size is chosen according to the above LHC integrated luminosity convention, which leaves the high-mass 
scenario m x = 1 TeV with essentially no background, and the low-mass scenario with about 24 background events in 
the m x mass window between 230 and 270 GeV [34], [35]]. The exact number of background events is not crucial for our 
analysis because we are mostly interested in understanding how well different signal hypotheses can be disentangled, 
rather than how well a signal can be separated from the background. 

We now imagine that 30 non-SM resonance-like events have been experimentally observed and we need to understand 
their origin. As the very first step, we may want to find out if the observed events are described by a resonance with 
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particular quantum numbers and a particular structure of interactions with the SM particles. An efficient way to do 
this is to compute a confidence level at which one hypothesis about resonance quantum numbers can be separated 
from another hypothesis. For example, we may ask if an observed resonance looks more like a SM Higgs boson or more 
like a graviton with minimal couplings to SM fields. We use our MC event simulation to check how well separation 
of the two hypotheses can be done in practice. Once we have a good idea about the nature of the resonance X, we 
can proceed to a more statistically demanding analysis whose goal is to determine all parameters which describe the 
angular distributions within a particular hypothesis. Then, one may attempt to interpret those parameters within a 
particular model of BSM physics. 

Note that the second approach clearly supersedes the first one and is the most optimal way to analyze data. Indeed, 
in the extreme version of the second approach, everything - mass, width, spin, parity, and polarization parameters 
of the particle X - is obtained through a single fitting process. The shortcoming of this approach is that it is very 
demanding statistically, since a large sample of events is required to enable simultaneous determination of a large 
number of signal parameters with decent precision. For the illustration of the second approach, we assume that 150 
signal events are available. 

We note that our methods are general enough to implement either of the two approaches described above. The 
only practical difference in application of these methods is that the polarization parameters are either fixed or left 
as free parameters in the fit. Since it is not practical to cover the continuous spectrum of possible scenarios in this 
paper, we choose several examples to illustrate these approaches. These scenarios are specified in Table HI 

Our studies are encouraged by the multivariate techniques developed for polarization studies in B-physics [4lj . 
where more than two dozen parameters describing angular and mass distributions are extracted in a simultaneous fit 
to a sample of several hundred events. We use an unbinned extended maximum-likelihood (ML) fit [U |42[ to extract 
simultaneously the signal and background yields, the background shape, and the parameters of the signal angular 
distributions. The likelihood function for N candidate events is written as 

C = cxp ^- Yynj - n hk ^j Y[ ( X] n J x Vj&i', Cj! £) + n bk s X Vbk s (xi; £)1 > ( 41 ) 

where n,, is the number of signal events for each resonance spin J, nbkg is the number of background events and 
V(xi] is the probability density function for signal or background. It is assumed that only one resonance is 
observed in a given mass window and the three yields, n,,, allow one to test different hypotheses. For this reason, no 
interference between the different resonances is considered. Alternatively, we can consider a signal event yield n s i g in 
place of the three n,, if we assume that mass distributions of resonances with different spins are identical and a single 
angular distribution can be written that incorporates angular distributions for all three spins as limiting cases. An 
example of the latter approach is given in Eq. (I38[) which describes the general spin-two distribution in four angles, 
but also includes both Eqs. (|35p and (|37p written for spin-zero and spin-one. 

Each event candidate i is characterized by a set of six observables Xi = {m zz , cos 6*, $i, cos6q, cos 62, The 
number of observables can be extended or reduced, depending on the desired fit. The signal polarization parameters 
{/aiA 2 ; < / , AiA 2 ; fzm} are collectively denoted by C 7 , an d the remaining parameters by When several decay channels 
must be combined, the joint likelihood in Eq. (1411) is the natural way to achieve this. In the current analysis, a 
combination of all lepton channels with electrons and muons in the ZZ decays is assumed, but it can be easily 
extended to include hadronic decays of Z bosons as well. 

For the signal, the probability density function Vj(xi \ Cj] £) is taken to be a product of a function of the angles and 
a function of the resonance mass. For the background, the probability density function is a product of uncorrelated 
functions for each observable. The assumption of small correlations among the discriminating variables, except for 
angles where relevant, can be validated in the selected data sample by evaluating correlation coefficients and is further 
tested with generated experiments. As long as the results of the generated experiments are unbiased, small correlations 
are acceptable. Since the parameterization of the background probability density function also depends on detector 
effects, similarly to the signal, experiment-specific treatment will, eventually, be required. It is an approximation of 
our approach that the background angular distributions are not correlated, which we find acceptable for the current 
analysis. 

The angular probability density function for the signal is the ideal distribution in Eq. (|28[) . multiplied by an 
empirically determined acceptance function <?(cos 6*, $1, cos#i, cos 62, Y x ), which can be evaluated independently 
for each event with the X particle rapidity Y x . We approximate the acceptance function by Q = Q^(cos 9*)x G^ 2 \&i), 
neglecting all correlations between different angles; the justification for this approximation was given in the previous 
section. The parameterization of the signal acceptance functions Cy*- 1,2 -* is compared with simulated distributions 
in Fig. [3j The acceptance function constructed along these lines can be included in the overall probability density 
function. 

It is important to point out that one can do a better job if all that is needed is to understand how well a particular 
hypothesis describes experimental data. Since we test a specific hypothesis, ratios of helicity amplitudes are fixed. 
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FIG. 4: Distribution of 21n(£i/£2) with the likelihood C evaluated for two models k = 1,2 and shown for 1000 generated 
experiments with the MC events generated according to model one (k — 1, open dots) and model two (k = 2, solid dots). Left 
plot: + vs. 0~; right plot: + vs. 2^. Effective signal hypothesis separation power S is 4.1 (left plot) and 2.8 (right plot). 



Hence, it is possible to determine the likelihood function numerically through the MC simulation of the full five- 
dimensional angular distribution, including the detector effects, though averaged over the Y x values. Effectively, this 
requires storing a normalized five-dimensional histogram with fine binning. Dependence on Y x could in principle be 
tabulated as well. A similar approach can be used for background parameterization. 

However, it is more difficult to employ MC computation of the signal probability density, if we want to treat all 
the parameters, including helicity fractions, as free parameters in the fit. A straightforward use of the MC program 
would require generating and storing probability distributions for a large set of independent polarization parameters. 
Clearly, for a large enough number of parameters to be fitted, this procedure becomes impractical. Instead, one could 
factorize the angular distribution into a sum of simple terms expressed as the product of a function of the angles and a 
function of the polarization parameters, as we have illustrated for the spin-zero case in Eq. (|35p . A single independent 
MC simulation for the angular dependencies of each such term can be performed and the results can be combined 
to produce a normalized angular parameterization with free parameters. However, for the purpose of this paper, we 
find the simplified approach to be acceptable since, as we have verified with our MC program, biases in the results 
are small. 

The only other kinematic observable that is used in the fit in addition to the angles is the ZZ invariant mass m zz . 
The background m zz distribution is parameterized with a polynomial and its shape can be left unconstrained in the 
fit. We make the following assumptions about the mass and the width of the X resonance. We assume that its width 
is negligible compared to the detector resolution effects, which are estimated to be about 3.5 GeV at m x = 250 GeV 
and 23 GeV at m x = 1 TeV with our simulation. The exact value of the resolution is only relevant for signal and 
background separation at the lower mass; the values quoted above are typical for both muon and electron final states. 
Hence, we describe the signal mass distribution with a Gaussian function of a known width. Should the resonance have 
a finite width that is larger than detector resolution, we can incorporate it into our analysis by replacing the Gaussian 
function with a convolution of the Breit-Wigner and Gaussian distributions. However, since the mass distribution has 
little, if any, influence on the resonance spin and its couplings to SM fields, a detailed study of mass resolution issues 
is outside the scope of this paper. 

The mass of the resonance m x is assumed to be known and is fixed in the fit. This assumption only affects the 
signal and background separation and requires certain care in the interpretation of the discovery significance. Detailed 
study of this interpretation is also outside the scope of this paper, but we note that, as Figj2]suggests, angular analysis 
improves the signal-background separation. As an example, consider 30 events due to the SM Higgs-like resonance at 
250 GeV. If only m zz is employed in the fit, the measured signal significance is 5.7 a and it increases by nearly twenty 
percent if also angular variables are included in the fit. Here significance is calculated from 21n(£i/£2) using the 
likelihood ratio for two scenarios with signal and with pure background, analogous to the technique discussed below, 
and corresponds to the Gaussian probability for pure background to fluctuate to the observed signal yield. The only 
constraint on background size comes from data sidebands included in the fit. 

We are now in position to discuss the results of our analysis. As an illustration of the first approach, in Fig. [H 
we plot the quantity 21n(£i/£2) for two sets of MC experiments generated with the models + and 0" (left plot) or 
+ and 2+ (right plot), as listed in Table [I] Each experiment has on average 24 background and 30 signal events, 
Poisson distributed. These events are generated at m x = 250 GeV in 1000 statistically independent experiments. 
The likelihood Ck in Eq. (|4Tj) is evaluated independently for each hypothesis k. The discussion of the statistical 
interpretation of 21n(£i/£2) was given recently in Ref. [l7|. The distribution of 21n(£i/£2) is expected to peak at 
positive values for events generated according to model one (on average C\ > £2) and at negative values for events 
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TABLE II: Results of the hypothesis separation of the seven scenarios listed in Table U for m x = 250 GeV. The five numbers 
quoted in each case correspond to the 1D/2D/3D/4D/5D angular information used in the fit, as discussed in the text. The 
results are reported as an equivalent separation of two single-width Gaussian distributions in terms of the number of Gaussian 
standard deviations between the two peaks. 



rr 1+ ir 2^ i\ sr 

0+ 0.0/0.0/3.9/4.1/4.1 0.8/1.0/1.8/1.9/2.3 0.9/1.0/2.5/2.6/2.6 0.8/0.9/2.4/2.5/2.8 2.6/2.6/0.0/2.6/2.6 1.6/1.7/2.4/3.0/3.3 
0" 0.8/1.2/2.8/3.0/3.1 0.9/1.0/2.5/2.8/3.0 0.8/0.8/1.7/2.0/2.4 2.9/2.9/4.1/4.8/4.8 1.6/1.7/2.0/2.7/2.9 
1+ 0.0/1.1/1.1/1.2/2.2 0.1/1.2/1.3/1.4/2.6 2.8/2.8/1.9/3.5/3.6 2.5/2.4/1.2/2.7/2.9 
1" 0.1/0.1/1.3/1.5/1.8 2.8/2.9/2.5/3.8/3.8 2.5/2.6/0.6/2.8/3.4 
2+ 2.9/2.9/2.6/3.6/3.8 2.3/2.5/0.5/2.5/3.2 
2j 3.6/3.6/2.5/4.2/4.3 



TABLE III: Results of the hypothesis separation of the seven scenarios for m x = 1 TeV with the same notation as in Table ITT1 

Q- 1 + 1" 2+ 2+ 2~ 

0+ 0.0/0.0/3.9/4.0/4.1 0.4/0.6/2.2/2.2/2.4 0.3/0.5/2.2/2.1/2.4 1.2/1.1/2.8/2.7/2.8 1.0/1.0/0.0/1.2/1.2 1.2/1.0/2.0/2.7/2.8 

0" 0.5/0.9/2.8/2.9/3.1 0.5/0.8/2.7/2.9/3.1 1.4/1.5/1.4/2.1/2.1 1.1/1.1/4.0/4.3/4.5 0.8/0.9/2.9/2.8/2.9 

1+ 0.0/1.1/1.2/1.2/2.3 0.9/0.9/2.0/2.0/2.3 1.3/1.4/2.1/2.4/2.6 1.3/1.3/1.1/1.8/2.0 

1" 0.9/0.9/2.0/2.0/2.1 1.3/1.3/2.2/2.5/2.7 1.4/1.6/0.0/1.4/2.2 

2+ 1.7/1.7/3.0/3.2/3.2 2.1/2.0/2.1/2.6/2.6 

2j 1.5/1.5/2.1/3.3/3.3 



generated according to model two (on average £2 > From Fig. [4j we extract the quantity S in the following way. 
We find the point beyond which the right-side tail of the left histogram and the left-side tail of the right histogram 
have equal areas. These areas correspond to the one-sided Gaussian probability outside of the S/2 a range. If the 
two histograms in Fig. [4] were perfectly Gaussian distributed with unit width, then S corresponds to the separation 
between the peaks of the two distributions. 

The above procedure is repeated for all combinations of seven hypotheses listed in Table UJ for m x — 250 GeV and 
m x = 1 TeV, and results are presented in Tables [TTI and Hill In those tables, we show the increase in the separation 
power between two hypotheses if more information about angular distributions is included in the likelihood fit. It 
is natural to refer to a likelihood fit that includes n angular variables as n-dimensional or nD; the results of the 
ID, 2D, 3D, 4D, and 5D fits are shown in Tables [TTI and IIIII In the ID fit, only the cos 0* angular distribution is 
used; in the 2D fit, the (cos0*,$i) distribution is included; in the 3D fit, we have (cos 01, cos02, $)i in the 4D fit, 
the (cos #*, cos 0i, cos #2,$) distribution is used; and finally in 5D fit, the (cos0*, $1, cos0i, cos0 2 , $) distribution is 
employed. Statistical significance achieved in the 4D fit is already close to the statistical significance of the optimal 
analysis which employs all five angles. We note that because we generate a finite number of experiments (1000), the 
statistical uncertainties on the values of S reported in Tables [TTI and [TTT1 are typically in the range 0.1 — 0.2. This 
explains some of the small discrepancies in the values quoted. For example, one could expect the separation of + 
and 0~ hypotheses to be the same in 3D, 4D, and 5D cases. However, statistical fluctuations and better background 
suppression make the 4D and 5D cases look better. 

It is clear from Tables [TT] and IIIII that the multidimensional angular analysis allows efficient separation of different 
spin hypotheses and, therefore, is a very powerful approach to spin determination. While separations 5 = 2 are 
typically achieved, there are important cases when larger separations are possible. In particular, we observe 5 = 4 
for the separation of + and 0~ hypotheses, which might be useful even for the SM Higgs boson if its mass exceeds 
2m z . Also, the separation of the spin-two with the minimal couplings hypothesis from the spin- two with longitudinal 
amplitude dominance hypothesis is close to S — 3, for any production mechanism scenario. In all cases, the results of 
the 5D fit provide the best separation power and supersede the other fits. Inclusion of the $1 angle in the 2D or 5D fit, 
compared to the ID or 4D fit, usually adds little information, except for the separation of the two spin-one hypotheses 
which have identical cos 0* distributions. Therefore, the separation power obtained from the angular analysis with all 
five angles is very close to the separation power of the 4D fit. 

A glance at Tables [TTJ and IIIII is sufficient to recognize that there are cases when different hypotheses about the 
nature of the resonance cannot be separated (5 = 0). For example, the production angles cos0* and $1 provide no 
separation between the two spin-zero scenarios (0 + and _ ) due to the absence of spin correlations. At the same 
time, the three helicity angle (cos 0i, cos 02, $) distributions are identical for the + and 2 \ hypotheses, therefore the 
3D case results in S = as well. The only source of separation in this case is the production angle. We have not 
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FIG. 5: Top: distribution of the number of fitted signal events n s i g (left) and the fraction of transverse component in the 
decay amplitude (/++ + /--) (right) in 1000 generated experiments with + hypothesis corresponding to Table ITVl Bottom: 
distribution of the above parameters normalized by the fit errors. 

considered any example with no separation from either production or helicity angles, but such situation is possible 
when spin-zero and spin-two scenarios with f z \ = f Z 2 = 0.4 are considered, as we discussed earlier. A hint at this 
situation is observed in the decreased 5-value that characterizes the separation of + and l\ hypotheses when the 
mass of the resonance increases from 250 GeV to 1 TeV. This is the consequence of the fact that at 1 TeV we allow 
qq production. This leads to j z \ — 0.25 and f Z 2 = 0.30, which is closer to the unpolarized case. As the result, the 
5-value decreases dramatically. 

Obviously, in the above studies we could not cover all possible scenarios. Just as an example, one could consider 
a scenario of a non-SM Higgs-like scalar with quantum numbers + with a\ = 0, a 2 7^ in Eq. ([2]). The separation 
of this hypothesis from the SM Higgs-like scalar will depend strongly on the mass of the resonance. As it is evident 
from Eq. (|16p . the single ^4oo amplitude dominates for m 2 x /rn^ ^> 1 in both cases. Therefore, separation of the two 

hypotheses is impossible at high mass. However, at lower mass, there is a sizable contribution of both A ++ and A 

in the SM Higgs-like case, but not in the other case. Therefore, the most optimal analysis strategy is to fit for all the 
polarization parameters, and then interpret the result within a particular model of BSM physics. We describe the 
feasibility of such fit in the remainder of this section. 

The most general analysis of an observed resonance requires that the likelihood in Eq. (|41[) is maximized with 
respect to the signal yield and all signal angular parameters, £ 7 , as well as with respect to unconstrained background 
parameters £. This requires a large number of free parameters even for signal alone. The fit delivers the most 
probable values of these parameters and the covariance matrix which describes their uncertainties and correlations 3 . 
In principle, it is possible to pursue this strategy within our approach, but for the fit to converge with high success 
rate, one may need large enough event sample. Therefore, for the illustration of the fit technique, we consider a 
limited procedure where only one spin hypothesis is considered at a time, but all other parameters relevant for that 
particular spin hypothesis are allowed to float in the fit. 

In Tables HVl and fVl this technique is illustrated with the spin-zero hypothesis with samples generated according to 
the + and 0~ scenarios, respectively. We show the results of the fit with both "perfect" and realistic detectors. In 
Fig. [5] we show distributions for two parameters in this analysis for a set of 1000 generated experiments. The means fj, 
and widths w of these distributions are shown in Table HVl as n±w. These values are the most probable central values 
and errors of a measurement, but they are subject to statistical fluctuation in each given experiment. Normalized 
distributions in Fig. [5] illustrate proper error estimates in the fit. In Tables [VTl and IVTfl the result of a similar study 
with the spin-one hypotheses 1 + and l - are presented. Finally, in Tables IVII1 IVIII1 and IIXI we illustrate the fit 



3 Knowledge of the full dependence of —2 In C on all the parameters is even more beneficial. 
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TABLE IV: Results of the fit for the free parameters of the spin-zero hypothesis with generated samples of SM Higgs-like X 
corresponding to + in Table Ull Experiments have been generated with two X masses according to two models in each case, 
with and without detector effects. 

mx = 250 GeV m x = 1 TeV 

generated fitted generated fitted 

without detector with detector without detector with detector 







150 


150 ± 13 


153 ± 15 


150 


150 ± 12 


152 ± 12 


(/++ 


+ /-) 


0.208 


0.21 ±0.07 


0.23 ± 0.08 


0.000 


0.00 ±0.03 


0.00 ± 0.03 


(/++ 


"/— ) 


0.000 


0.01 ±0.13 


0.01 ±0.14 


0.000 


0.00 ± 0.02 


0.00 ± 0.02 


(0++ 


±0— ) 


2tt 


6.30 ± 1.46 


6.39 ± 1.54 


2tt 


free 


free 


(0+ + 


-0— ) 





0.00 ± 1.06 


0.01 ± 1.09 





free 


free 



TABLE V: Results of the fit for the free parameters of the spin-zero hypothesis with generated samples of a pseudo-scalar X 
corresponding to 0~ as discussed in Table [TV] 







generated 


m x = 250 GeV 
fitted 

without detector with detector 


generated 


m x = 1 TeV 
fitted 

without detector with detector 






150 


150 ± 13 


151 ± 15 


150 


151 ± 12 


150 ± 13 


(/++ 


+ /— ) 


1.000 


1.00 ±0.05 


1.00 ±0.06 


1.000 


1.00 ±0.05 


1.00 ± 0.06 


(/++ 


-/-) 


0.000 


0.00 ±0.35 


0.00 ± 0.40 


0.000 


0.00 ±0.31 


-0.01 ±0.32 


(0+ + 


±0— ) 


N/A 


free 


free 


N/A 


free 


free 


(0+ + 




7T 


3.15 ±0.31 


3.14 ±0.41 


7T 


3.15 ±0.31 


3.14 ±0.33 



TABLE VI: Results of the fit for the free parameters of the spin-one hypothesis with generated samples of an exotic pseudo- 
vector X corresponding to 1 + as discussed in Table ItVl 







m x = 250 GeV 




m x = 1 TeV 




generated 


fitted 


generated 


fitted 






without detector with detector 




without detector with detector 




150 


150 ± 13 152 ± 15 


150 


151 ±12 152 ±13 


/+o 


0.250 


0.26 ± 0.15 0.26 ±0.17 


0.250 


0.25 ±0.16 0.25 ±0.15 


O+o - 


0o-) 


0.02 ± 0.73 0.01 ±0.86 





-0.03 ±0.67 -0.03 ±0.77 



TABLE VII: Results of the fit for the free parameters of the spin-one hypothesis with generated samples of an exotic vector 
X corresponding to 1~ as discussed in Table ITVl 







mx = 250 GeV 




m x = 1 TeV 




generated 


fitted 


generated 


fitted 






without detector with detector 




without detector with detector 


?l s ig 


150 


150 ±13 152 ±16 


150 


151 ±12 152 ±13 


/+0 


0.250 


0.24 ± 0.15 0.24 ±0.16 


0.250 


0.24 ±0.17 0.25 ±0.15 


(0+0 - 


00-) 7T 


3.14 ± 0.72 3.18 ±0.76 


7T 


3.14 ±0.71 3.14 ±0.69 



with the spin- two hypotheses for 2+, 2 J, and 2 _ , respectively, where we also restrict the hypothesis to the definite 
parity to reduce the number of free parameters. We used 1000 generated samples of 150 signal events each, Poisson 
distributed, and the corresponding number of background events, per sample. In all cases the set of observables 
includes m zz , cos 9*, cos 0±, cos 82, and There are several cases where contribution of a certain helicity amplitude 
is either exactly zero (e.g. Aqo in the 0~ case) or negligible (e.g. ^4++ in the 2+ case). Since measuring the phase of 
that amplitude becomes impossible, we leave it as "free" in the tables. 

Coming back to the question of separating SM-like and non-SM-like + hypotheses, it is clear from Fig. [5] and 
Table ITVl that the hypothesis 04 = and a-2 ^ in Eq. ([2]), corresponding to /++ ± /__ = 0, can be separated 
with high confidence at 250 GeV, but no separation is possible at 1 TeV due to the absence of a sizable transverse 
amplitude in both scenarios. As another example, the separation of the 1 + and 1~ scenarios in Tables HT1 and ITTTl is 
related to the measurement of (0+o — 4>o~) in Tables IVTl and rVTfl although there is an additional free parameter in the 
latter case. The signature of the minimal coupling model of a spin-two resonance decay is the dominance of the 
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TABLE VIII: Results of the fit for the free parameters of the spin-two hypothesis with positive parity with generated samples 
of an exotic X corresponding to 2+ as discussed in Table ITVl 





generated 


m x = 250 GeV 
fitted 

without detector with detector 


generated 


m x =l TeV 
fitted 

without detector with detector 


n s ig 


loU 


150 ± 13 


151 ± 16 


loU 


151 ± 12 


153 ± 13 


fz2 


1.000 


1.00 ±0.17 


0.84 ±0.17 


0.750 


0.75 ±0.12 


0.80 ±0.10 


hi 


0.000 


0.00 ±0.19 


0.00 ±0.25 


0.250 


0.25 ±0.14 


0.16 ±0.15 


/++ 


0.013 


0.01 ±0.04 


0.00 ±0.05 


0.000 


0.00 ±0.05 


0.00 ± 0.05 


/+- 


0.282 


0.28 ±0.04 


0.31 ±0.05 


0.445 


0.44 ± 0.06 


0.44 ± 0.04 


f+o 


0.075 


0.07 ± 0.04 


0.06 ± 0.05 


0.000 


0.01 ±0.06 


0.01 ±0.06 


<t>++ 





0.00 ± 1.75 


0.04 ± 1.76 





free 


free 



TABLE IX: Results of the fit for the free parameters of the spin-two hypothesis with positive parity with generated samples 
of an exotic X corresponding to 2 J as discussed in Table IIVI 







m x = 250 GeV 






m x = 1 TeV 






generated 


fitted 




generated 


fitted 








without detector with detector 




without detector with detector 




150 


150 ± 13 


154 ± 15 


150 


151 ± 12 


152 ± 13 


fz2 


0.400 


0.40 ±0.07 


0.33 ±0.10 


0.300 


0.30 ±0.08 


0.34 ± 0.09 


f*l 


0.000 


0.00 ± 0.03 


0.00 ± 0.04 


0.250 


0.25 ±0.07 


0.23 ± 0.07 


f++ 


0.104 


0.10 ±0.04 


0.06 ± 0.05 


0.000 


0.00 ±0.01 


0.00 ± 0.02 


f+- 


0.000 


0.00 ± 0.04 


0.04 ± 0.05 


0.000 


0.00 ±0.01 


0.00 ± 0.03 


f+o 


0.000 


0.00 ± 0.02 


0.00 ±0.02 


0.000 


0.00 ±0.02 


0.00 ± 0.02 


4>++ 


7T 


3.20 ±0.75 


3.17 ±0.71 


7T 


free 


free 



TABLE X: Results of the fit for the free parameters of the spin-two hypothesis with negative parity with generated samples 
of an exotic X corresponding to 2 _ as discussed in Table ITVl 





generated 


m x = 250 GeV 
fitted 

without detector with detector 


generated 


m x = 1 TeV 
fitted 

without detector with detector 


71 s ig 


150 


150 ± 13 


152 ± 15 


150 


151 ± 12 


152 ± 13 




1.000 


1.00 ± 0.06 


0.97 ±0.12 


0.750 


0.75 ±0.11 


0.81 ±0.15 




0.000 


0.00 ±0.11 


0.05 ±0.20 


0.250 


0.25 ±0.13 


0.19 ±0.24 


/++ 


0.125 


0.12 ±0.06 


0.13 ±0.07 


0.000 


0.00 ± 0.07 


0.00 ± 0.07 



and A (- amplitudes, while the model discussed in Refs. [2ll - l23| predicts longitudinal Aoo amplitude dominance. This 

was our motivation for comparing the 2+ and 2 \ cases. We observe that the angular distributions of the two models 

could be well separated based on the measurement of /-| , as it is also evident from Fig. [2] and Tables HH HID There 

are several cases where probability distributions of some parameter are not Gaussian; typically, this occurs when the 
central value of a parameter is at, or very close to, its physical boundary. In those cases, we still approximate the 
expectation value for this parameter and the error by y, ± w, where we extrapolate the value of w from the rms, but 
we note that statistical meaning of the error estimated in such a way can be very different. For this reason, the proper 
interpretation of the values quoted in the tables in terms of the separation S value requires a more careful treatment 
which accounts for the non-Gaussian shape of all distributions. 

The general conclusion from Tables IIVI [Xl is that our analysis yields unbiased results for all the parameters that 
describe the angular distributions in the production and decay of a resonance X. Indeed, it follows from those tables 
that, when detector effects are neglected, the central values of all the parameters correspond to their input values. 
This verifies both the fit implementation and angular parameterizations. 

Once detector effects are introduced, certain biases appear but they are minor. These biases come from the fact 
that we ignore angular correlations in both the detector acceptance function and in our parametrization of angular 
background distributions. Both of these parameterizations are detector-specific and would need to be modeled in each 
experiment. For illustration purposes, we introduced a simplified model of the acceptance function in the present 
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paper. We pointed out that the sophistication of the required model for the acceptance function and the background 
parametrization depends on the actual mass of the observed resonance. Indeed, it is evident from Tables ITVHXl that the 
bias is smaller at higher masses due to reduced angular correlations, as we explained earlier. The expected statistical 
errors of the measurements are not strongly affected by detector effects; this allows us to predict statistical precision 
in future measurements. 

The above studies illustrate the precision on individual parameters that can be achieved with a sample of 150 signal 
events. These results can be easily extrapolated to an event sample of a different size. On the other hand, if an 
accurate precision estimate is required in a particular scenario, dedicated analysis is necessary. In the present paper, 
we do not include the evaluation of systematic uncertainties beyond simplified studies of resolution and acceptance 
detector effects, but these and other systematic effects are typically much smaller than statistical uncertainties if the 
signal significance is close to the discovery threshold. 

VII. CONCLUSIONS 

The importance of the LHC physics program for the future of high-energy physics requires designing techniques 
and analysis tools that are not biased towards a particular model of beyond the Standard Model physics. While it 
may be impractical or impossible to do so in many complicated cases, the case of a single resonance production at 
the LHC is sufficiently robust to aim at the most general description, yet simple enough to make such a description 
possible. Motivated by this, we studied production of a single resonance at the LHC and its decay into four-lepton 
final states in the process pp — > X — » ZZ — » hhhh, allowing for the most general couplings of the resonance to 
Standard Model matter and gauge fields. Our goal was to understand if full reconstruction of the final state and the 
analysis of the most general angular distributions can be used to determine the spin and parity of the resonance and 
to constrain its couplings to vector gauge bosons, quarks, and gluons. 

It is interesting to point out that, in some cases, generality of the couplings leads to new effects in the angular 
distributions. For example, imagine that the spin-two resonance is identified with the massive graviton whose coupling 
to gauge and matter fields is fixed through the energy-momentum tensor. Production of such a resonance in gluon 
fusion occurs only with spin projections on the collision axis J z = ±2. On the other hand, if a general coupling of the 
spin-two resonance to gluons is considered, production with J z = becomes possible as well. There are interesting 
consequences of this observation related to a possibility to discriminate between minimal and non-minimal couplings 
of the spin-two resonance through the analysis of angular distributions. 

Motivated by this example, we derived the most general angular distributions of four fermions in the process 
pp — > X — > VV — > /1/1/2/2, considering spin-zero, spin-one, and spin-two resonances . We used those angular 
distributions to construct the likelihood functions for the angular analysis. We wrote a Monte Carlo simulation 
program that describes the production of the resonance X in the process pp — > X — > ZZ — > hhhh, includes all spin 
correlations throughout the decay chains, and employs the most general couplings of the resonance X to matter and 
gauge fields of the Standard Model. We supplemented this analysis with a simplified (but fairly realistic) model of 
the detector effects. When all the pieces of our study are put together, we obtain a powerful analysis tool that allows 
for a realistic estimate of how much information the two LHC experiments will be able to extract from the study of 
angular distributions once the resonance in the process pp — > X — > ZZ hhhh is observed. 

We find that angular distributions provide good separation of various hypotheses about the spin of the resonance. 
As we illustrate in Tables |H] and IIII1 with 30 fully reconstructed events a typical separation of various spin hypotheses 
is S = 2 — 3 er, as defined in text, but in some cases separation as good as 4 a is achieved. With a somewhat larger 
event sample, it becomes possible to determine helicity amplitude fractions that characterize resonance production and 
decay as well as the resonance spin from a multidimensional likelihood fit with decent precision. Model-independent 
determination of the helicity amplitude fractions and phases is the ultimate goal of such an analysis which can then be 
interpreted within any model of beyond the Standard Model physics. We provide relationships between the amplitude 
measurements and the fundamental coupling constants of the resonance to matter and gauge fields. Our studies show 
that such model-independent analysis is viable at the LHC. We look forward to its application to real LHC data. 

Finally, we note that the analysis presented in this paper can be extended in a number of ways. A natural possibility 
is to allow hadronic final states (e.g. Z — > qq) and/or missing energy. An obvious candidate for the latter is the 
decay X — > W + W~ . We note that most of the discussion of the X — > ZZ decays given in this paper also applies to 
decays X —> W + W~ . However, there are important differences related to the fact that W + and W~ are not identical 
particles. Consequently, Eq. ([14")) does not apply in the W + W~ case and nine, rather than six, independent helicity 



4 Angular distributions for two-particle decays of the X bosons are given in the appendix. 
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amplitudes are required to describe X — > W + W~ decay. While the increase in the number of the helicity amplitudes 
implies that most of the formulas for angular distributions derived in this paper are not complete for X — > W + W~, 
it is interesting to remark that Eq. ([2]) remains the most general description of a spin-zero X — > W + W~ decay and 
Eq. (J3SJ) is valid, with Ri = R 2 = 1. 

Another aspect of the decay of the resonance X to two non-identical particles is that forward-backward asymmetry 
can be generated. Of course, this requires that the initial state is asymmetric (e.g. qq) as well. In general, the 
forward-backward asymmetry manifests itself in the odd terms in cos 9* distribution. To isolate and measure these 
terms, an unambiguous definition of the direction of the z-axis is required; a suggestion on how to do that was given 
in Ref. [43j]. On the experimental side, the angular analysis of X — ► W + W~ — > 41 decay is challenging due to the 
presence of two neutrinos in the final state. It might be beneficial to consider semileptonic final states where one W 
decays hadronically and the other leptonically. However, this case requires detailed studies because of potentially large 
irreducible backgrounds from events with jets and missing energy. Similar issues should be investigated if semileptonic 
final states are allowed in X — » ZZ decays. We hope to return to the discussion of these channels in the near future. 
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Appendix A: Angular distributions for other decay channels 

In this appendix, the angular distributions for a resonance decaying into two particles, including 77, gg, l + l~ , and 
tt, are given. The general formula for the two-particle production and decay gg or qq — > X — > P\P 2 can be obtained 
from Eq. ([25| by changing the spin-quantization axis of the X to z and setting fl = ($1, 9* , — $1) and ft* = (0, 0, 0), 
where $1 is now an arbitrary azimuthal angle to be integrated out. The resulting formula reads 



dT(Xj -► P 1 P 2 ) = / _ 1 



rdcosfl* V 2 



E a^E/'^U-a^*)) 2 , (ai) 



Ai .A2 



where Ai and A2 run over all possible helicities of P\ and P 2l m runs over all possible X spin projections, and 
/aiA 2 = I^AiA 2 | 2 /X) l^fcz| 2 ! as defined earlier for the X — » ZZ decay. However, this definition is more general and 
includes Ai = ±1/2 for decays to fermions X — > ff. In the latter case, we choose the same notation as in Eq. (|29|) 

for their parameterization with /++, / , /_| , / \-, where we omit the 1/2 for simplicity. Equation (IAI[) includes f m 

parameters, where /o = f z o, f±2 = /z2/2, and f±i = (f z i ± A/ z i)/2, allowing for a quark direction measurement in 
qq production and, generally, for two non-identical particles Pi and P 2 . 

We now specialize to particular channels. For the X — > 77 decay channel, the angular distribution can be obtained 

cither from Eq. (jAII) or by integrating out the $1 dependence in Eq. (|33|) . Note that only /++, / , and /_+ = /_| 

parameters are non-zero in the decay to two massless photons. Therefore 

16dr(A,^ 2 ->77) = (2 _ 2U + fz2) _ 6(2 _ ^ _ fz2) cog2 r + 3(6 _ _ 5/zz) cog 4 r 

51a cos 0* 

+/+_ {(2 + 2/, 1 -7/ 22 ) + 6(2-6/ zl + A 2 )cos 2 r-5(6-10/ 2l -5/ 22 )cos 4 r} . (A2) 

This equation describes the most general case. The special case of the minimal coupling in both production and decay 

corresponds to f z \ + f z2 = 1 and = /_+ = 1/2. In this case, one obtains (I + 6cos 2 #* + cos 4 0*) for the gg 

production mechanism with f z2 = 1 and (1 — cos 4 #*) for the qq production mechanism with f z \ = I. Equation ()A2|) 
is also applicable to the decay to two gluon jets Xj =2 — > gg, but additional constraints could be used when the 

production and decay mechanisms are the same: (I — f z \ — f z2 )/f z2 = (1 — f-\ — /-+)/(/h + / h)- 

For the decay to a fermion-antifermion pair, we obtain 

S ' /J ( A ' " ' 7 = (/+- + /-+)(! + cos 2 9*) + 2 (1 - /+_ - /_+)(! - cos 2 9*) + 2 (/+_ - f-+)Af zl cos9* , (A3) 



3Td cos 9* 
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16dr 5Td^e* fI) = 2 {f+ - + f - +) { {hl + Iz2) + 3(2 ~ 3/21 " 2fz2) cos2 6 * ~ (6 " 10/21 " 5fz2) cos4 r } 

+ (1 - /+- - /-+) { (2 - 2f zl + f z2 ) - 6(2 - 4/,i - /, 2 ) cos 2 0* + 3(6 - 10f zl - 5/ 22 ) cos 4 0*} 
-4(/+_ - /_+)A/,i(cos0* - 2 cos 3 r) . (A4) 

where for a massless fermion in the final state (/++ + /--) = (1 — /+- — /-+) = 0, which would describe the decay 
X —> . It follows from this formula that there is a forward-backward asymmetry in this decay, as was pointed 
out in Ref. [l3[ in the context of spin-one decays to a fermion pair. A dilution factor needs to be introduced in front 
of the A/zi terms, which depends on the ability to measure the sign of cos 0* in an experiment. The special case of 
the minimal coupling in gluon fusion corresponds to f z ± + f Z 2 = 1. 



Appendix B: Supporting material 



Supporting material for this analysis may be found in Ref. [44j . where we provide the Monte Carlo simulation 
program and the most general angular distributions used in this analysis. For completeness, we present the general 
angular distribution in the production and decay of a spin-J particle X in parton collisions ab — > X — > ZZ — > 
(/l/iX/a/a)- In order to simplify expressions, we redefine the fifth angle from f 1 to $ = $1 + 3>/2, which can be 
interpreted as the angle between the production plane and the average between the two decay planes shown in Fig. [T] 

AO dTj 



Yd cos 9*d^d cos 9\d cos 9 2 d$ 

F oo( 9 *) x { 4 /oo sin 2 9 X sin 2 9 2 + (/++ + /__) ((1 + cos 2 0i)(l + cos 2 9 2 ) + 4R X R 2 cos9 1 cos0 2 ) 



2(/-+ 



/__) (i?i cos 9i (1 + cos 2 



R 2 (l + cos 2 6»i ) cos ( 



+ 4y/f ++ f m (Ri - cos6>i) sin6»i(i? 2 - cos(9 2 ) s'm9 2 cos($ + <p ++ ) 

+ 4\/7 foo (Ri + cos 61 ) sin 6\ {R 2 + cos0 2 ) sin0 2 cos($ - <j> ) 

+ 2 sin 2 9 X sin 2 9 2 cos(2$ + 0++ - 

+4^(0*) x {(/ +0 + / _)(1 - cos 2 9 l cos 2 2 ) - (/+o - fo-)(Ri cos9 1 sin 2 9 2 + R 2 sin 2 0i cos0 2 ) 

+ 2 V/+o/o- sin 6»i sin6> 2 (i?i_R 2 - cos0i cos0 2 ) cos($ + cj> +0 - O -)j 
+(-l) J x 4Pf lt (d*) x {(/+„ +f _)(R 1 R 2 + cos0icos0 2 ) - (/ +0 - / o _)(i2icos0a + Jfccosfli) 

+ 2\//+o/o- sin 0i sin0 2 cos($ + </>+o — 0o-)| sin0i sin0 2 cos(2\&) 
+2F 2 i(0*) x + cos 2 0i)(l + cos 2 2 ) - AR X R 2 cos0i cos0 2 | 

+(-l) J x 2F J _ 22 (9*) x /+_ sin 2 9 1 sin 2 2 cos(4*) 
+2F J 2 (9*) x {2v^A7 e 



_ sin 1 sin W 2 x 



- COS 0i)(i?2 + COS 02 ) cos(2W - 



+ (Rl + COS0i)(i? 2 - cos0 2 )cos(2* + 
+ Vf++f+- sin2 'M 1 - 2i? 2 cos 2 + cos 2 2 ) cos(2* - $ 
+ (1 - 2Ri cos 0i + cos 2 0i ) sin 2 2 cos(2* + $ - 



+ sin 2 0i (1 + 2i? 2 cos 2 + cos 2 2 ) cos(2* - $ 

+ (1 + 2i?i cos 0i + cos 2 0i ) sin 2 2 cos (2* + $ + + _ 



>++) 
</>+- +0- 
</>-)]} 



-2V2F o i(0*) x {2v^/W^ 



sin0i(i?i - cos0i) sin 2 2 cos(* - $/2- < 
- sin 2 0i sin0 2 (i? 2 - cos0 2 )cos(* + $/2 + </> +0 ) 
+ 2^/oo/o- [sin 2 0i sin0 2 (i? 2 + cos0 2 ) cos(* + $/2 - <fo-) 
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- sin6»i(i?i + cos0i)sin 2 2 cos(* - $/2 + O _ 
+ VTf+7+o ( 1 - 2i?i cos 9 1 + cos 2 0i ) sin 2 (R 2 — cos 9 2 ) cos^ + <£>/2 + + _|_ — 4>+o) 

- sin 0i (.Ri - cos0i)(l - 2R 2 cos6 2 + cos 2 2 ) cos(* - $/2 - ++ + +o ) 
iin0i(i?i -cos0i)sin 2 2 cos(* + 3$/2 + ++ - <f> -) 



- sin 2 6»i sin6» 2 (i? 2 - cos0 2 )cos(* - 3$/2 - ++ + O _ , 

a 2 6 1 sin6» 2 (i? 2 + cos0 2 )cos(* - 3$/2 + 0__ - +o ) 



-sin0i(i?i + cos 6»i ) sin 2 2 cos(* + 3$/2 - + 0+o) 

+ y/f—fo- sin0i(it-i + cos0i)(l + 2i? 2 COS0 2 + COS 02 ) cos(^ — $/2 + 0__ — 00-) 

- (1 + 2R 1 cos 0i + cos 2 0i ) sin 2 (R 2 + cos 2 ) cos(* + $/2 - 0_ _ + O _ )] | 

-2V2 F{ 2 {6*) x {V/+-/+0 [(1 - 2i?i cos0i + cos 2 0i) sin0 2 (i? 2 + cos0 2 ) cos(* + $/2 - 0+_ + +o ) 

- sin0i(i?i +cos0i)(l - 2i? 2 cos 2 + cos 2 2 )cos(* - $/2 + 0+_ - +o ) 
+ y/f+-fo- sin0i(fli - cos0i)(l + 2i? 2 cos0 2 + cos 2 2 )cos(* - $/2 - 0+_ + O _) 



-(-1) J x 2V2 Fij 



sin0i(i?i - cos0i)(l + 2i? 2 cos 2 

- (1 + 2F1COS01 +cos 2 0i)sin0 2 (i? 2 - cos0 2 )cos(* + $/2 + 4>+- - 0o-)] } 
^12(0*) x { V/+-/+0 [sin0i(i?i - cos0i) sin 2 2 cos(3* + $/2 - 0+_ + +o ) 

- sin 2 0i sin0 2 (i? 2 - cos 2 ) cos(3* - $/2 + 0+_ - +o ) 



— sin 0i 



!-0+_ +0o_) 



- sin 2 0i sin 2 (i? 2 - cos 2 ) cos(3* - $/2 + 0+ 
;in 2 0i sin0 2 (i? 2 + cos 2 ) cos(3* - $/2 
?i + cos 0i ) sin 2 2 cos(3* + $/2 + + _ - O _) | , 

where Afj is the normalization factor which does not affect the angular distributions and the functions F-j(0*) 
defined as follows 

r;'(n= E f m d j mi (e*)d j mj (6*), (B2) 

m=0,±l,±2 

:, and f = f z0 = (1 - f zl - f z2 ). Note that for odd J one has / 00 = /++ = / — = 0, 
do not contribute. 



(Bl) 
are 



where f ±1 = f z i/2, f± 2 = /* 2 /2, and f = f z0 = ( 
and therefore F Q J J =odd (9*) terms do not contribute. 

Below we show F/j(6*) explicitly for J = 0, 1, and 2. For spin-zero, we have 

^oo = 1 . 

^11 - f-ll - ^22 — ^-22 — ^02 — ^01 - P 12 - *-12 — U ' 

where only four parameters/ ++ , /__, 0++, and 0__ remain relevant, and foo can be expressed as foo = (1 — /++ — 
/--). For spin-one, we have 



(B3) 



Fi 1 i = i (l + cos 2 0*), 
^11 = J(l-cos 2 0*), 

^oo — -^22 - ^-22 - -^02 - ^01 — ^12 — ^-12 — u ! 

ipressed as / _ = (1 - 2/ +0 )/2. 



UU -i^r Zl^r 

where only two parameters / +0 and [0 +o — 0o-] remain relevant, and / _ can be ex; 
Finally, for spin-two, we have 



(B4) 



Foo = l{( 2 ~ 2 + ~ 6(2 - 4/ z i - f z2 ) cos 2 0* + 3(6 - 10/,i - 5/ z2 ) cos 4 0*} , 
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Ffi = \{ihi + fz2) + 3(2 - 3/ 2l - 2/ z2 ) cos 2 0* - (6 - 10/ zl - 5/ z2 ) cos 4 0* } , 
J^U = ~\{(f*i ~ /»a) + ( 6 - 10 /,i - 5/-2) cos2 ^ } sin2 r . 

*aa = ^{(6 - 2U - 5 fa) - 6(2 - 2f zl - 3/ 22 ) cos 2 9* + (6 - 10/ 2l - 5/ z2 ) cos 4 6* } , 
^ 2 22 = ^{6 - 10/ 2l - 5/ z2 } sin 4 0* , 

^o 2 = " 2/»i ~ 3/ 22 ) - (6 - 10/,i - 5/, a ) cos 2 } sin 2 0* , 

*oi = -^{(2 - 4/,i - /, a ) - (6 - 10/ 2 i - 5/, a ) cos 2 0* } cosfl* sintf* , 
^12 = ^{(6 - 6/ 2i - 9/ z2 ) - (6 - 10/,! - 5/ z2 ) cos 2 0* } cosfl* sintf* , 

^-12 = ~ J (6 - 10/ z i - 5/ 22 ) cos 0* sin 3 0* , (B5) 
where all parameters contribute and again /oo can be expressed as /oo = (1 — /++ — / — %f+o — 2/o- — 2/_| ). 
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